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ABSTRACT 

The  presence  of  non— linear  behavior  in  potentiodynamic 
polarization  plots  has  resulted  in  difficulty  in  determinxnq 
the  Tafel  constants  from  such  plots.  A  FORTRAN  based  program 
involving  numerical  differentiation  techniques  using  a 
graphical  display  was  used  to  determine  the  existence  of  the 
Tafel  regions. 

Various  alloys  polarized  in  synthetic  seawater  and  a 
3.5/C  NaCl  solution  were  analyzed.   Although   severe 
concentration  polarization  often  dominated  the  cathodic 
branches  the  techniques  employed  did  allow  for  the  selection 
of  regions  which  approached  linear  behavior.   The  effects  of 
concentration  polarization  in  hindering  the  determination  of 
Tafel  constants  were  exemplified  by  the  uncovering  of  a 
cathodic  branch  containing  a  small  region  where  only 
actvivation  polarization  dominated  followed  by  the  onset  and 
total  domination  of  concentration  polarization - 

A  method  of  determining  where  the  anodic  and  cathodic 
currents  begin  to  dominate  the  potentiodynamic  polarization 
curve  is  introduced. 
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I •   INIRODUCJigN 

Recent  studies  at  the  Naval  Postgraduate  School  have 
centered  on  the  vibrational  characteristics  of  certain  high- 
damping  alloys,  [Refs.  1-8].   Most  of  these  studies  have 
been  concerned  with  either  the  measurement  of  the  damping 
capacities  of  these  alloys  and/or  with  the  heat  treatments 
necessary  to  attain  the  micros true tures  associated  with  the 
observed  high— damping  capacities-   In  addition,  two  separate 
studies  were  conducted  in  an  attempt  to  determine  the 
corrosion  characteristics  of  these  high  damping  alloys  in 
marine  environments. 

In  the  studies  involving  corrosion  characteristics,  both 
laboratory  and  actual  marine  environment  exposure 
experiments  were  conducted.   The  actual  marine  environment 
exposure  tests  calculated  corrosion  rates  based  on  direct 
weight  loss  methods.   The  laboratory  tests  conducted  were 
Linear  Polarization  and  Potentiodynamic  Polarization  in 
synthetic  seawater  [Ref.  *?]  and  a  3.5X  NaCl  solution  FRef. 
10]  environments.   Although  the  correlation  between  the 
laboratory  and  direct  exposure  methods  was  good,  the 
individuals  involved  in  both  of  these  laboratory  studies 
expressed  concern  with  the  difficulty  in  determining  the 
Tafel  constants  necessary  to  determine  corrosion  rates. 


The  actual  experiments  were  conducted  on  the  EG3tG 
Princeton  Applied  Research  Model  351  corrosion  measurement 
system.   The  Model  351  utilizes  a  PAR  Model  lOOO 
microcomputer  featuring  two  Motorola  68000  microprocessors 
and  a  touch— screen  input  to  create  an  environment  which  is 
well  suited  for  conducting  various  corrosion  tests.   The 
system  is,  however,  limited  by  the  lack  of  a  key  board  and 
more  importantly  by  an  operating  system  which  is  proprietary 
in  nature.   These  two  combine  to  make  any  analysis  of 
experimental  data,  other  than  by  the  means  provided  with  the 
operating  software,  impossible  unless  the  data  can  be 
transfered  to  another  computer. 


1 1 •   DATA  TRANSFER 

The  PAR  351  and  its  associated  software  has  the  ability 
both  to  plot  the  results  of  an    experiment  and  to  send  just 
the  data  points  to  a  printer  via  an    RS  232  serial  port. 
Although  the  software  supports  only  Houston  Instruments 
plotters,  the  printer  option  is  configurable  for  numerous 
printers  as  long  as  the  printer  supports  DCl  (XON,  ASCII  17) 
and  DC3  (XOFF,  ASCII  19)  protocol.   Other  RS  232 
requirements  such  as  baud  rate,  parity,  word  length,  and 
stop  bits  are    configurable  as  the  user  sees  fit. 

During  the  initial  data  transfers,  in  the  present  work, 
Hayes  Microcomputer  Inc.   SMARTCOM  II  modem  software  was 
used  to  capture  the  data.   This  method  was  successful  most 
of  the  time  but  since  the  two  systems  were  not  using  modems 
but  instead  were  directly  connected  via  their  serial  ports 
some  problems  were  experienced  and  editing  of  the  data  file 
after  transfer  was  always  necessary.   The  most  serious 
detraction  was  the  need  to  wade  through  various  menus  in 
order  to  accomplish  the  transfer. 

In  view  of  the  both  the  problems  experienced  with  using 
SMARTCOM  II  and  the  neccessity  for  the  user  to  be  familar 
with  the  way  SMARTCOM  II  works,  it  was  decided  to  write  a 


data  transfer/capture  program  which  would  both  speed 
execution  and  eliminate  the  need  for  user  knowledge  of 
communications  software. 

Using  the  sample  TTY  data  transfer  program  included  in 
the  manual  for  Microsoft  GW-BASIC  version  2  CRef.  11]  as  a 
core  a  BASIC  program  was  written  to  accomplish  the  transfer. 
The  core  was  modified  to  eliminate  the  transfer  portion  and 
most  importantly  line  feeds  and  carriage  returns  sent  by  the 
351  were  automatically  deleted  and  a  carriage  return  was 
inserted  only  at  the  end  of  a  line  containing  data.   This 
eliminated  the  problem  of  blank  lines  appearing  within  the 
data  file  corresponding  to  page  headers  and  footers.   Other 
modifications  which  were  made  involved  checking  for  and 
correcting  two  formatting  variations  which  resulted  in 
compatibility  problems  between  the  transfered  file  and  its 
use  in  FORTRAN  application  programs.   The  first  involved  a 
change  from  a  floating  format  to  a  exponential  format  when 
the  potential  was  equal  to  O.OOOO  volts.   The  351  normally 
sends  potential  with  4  significant  digits.   When  the 
potential  is  O.OOOO  volts  it  is  displayed  as  -XXE-12.   The 
formatting  shift  is  not  compatible  with  any  FORTRAN  formats 
as  the  total  field  width  changes  in  addition  to  the  minimum 
number  of  available  digits.   This  occurance  is  now  trapped 
and  the  file  modified  to  reflect  the  actual  potential  of 
O.OOOO  volts.   The  second  formatting  variation  always 
occured  and  did  so  in  a  region  which  was  of  extreme  interest 


lO 


as  follows.   When  the  current  density  changes  from  negative 
to  positive  there  is  a  range  where  the  exponent  is  E— 12. 
Unlike  FORTRAN  write  statements  which  always  reserve  a  space 
for  the  +  or  -  sign,  ie.  the  +  sign  is  not  printed  unless 
specified,  the  351  used  the  space  for  the  -  sign  if  present 
or  for  the  first  digit  if  not.   Again  this  resulted  in  a 
change  in  total  field  width  at  the  occurance  of  -a.bcdE— 12. 
The  FORTRAN  read  statment  would  read  this  to  be  -a.bcdE-1. 
This  was  handled  by  simply  deleting  the  least  significant 
digit  such  that  the  number  is  read  as  — a.bcE— 12  when  this 
occurs. 

In  its  final  form  the  data  transfer  program  requires  the 
user  to  reply  only  a  to  a  (R)eceive  or  (E)xit  prompt  and  to 
enter  a  file  name  if  the  (R)eceive  option  is  selected. 
Because  the  351  can  not  write  a  selected  configuration  to 
disk  the  user  is  still  required  to  properly  configure  the 
351  to  the  proper  settings  when  powering  up.   These  ^rez 


PROTOCOL      -  PRINT 
PARITY        -  NONE 
BAUD  RATE    -  2400 
STOP  BIT     -  1 
WORD  LENGTH  -  8 


Once  these  options  ^re    properly  set  on  the  351  the 
desired  data  file  should  be  loaded  into  the  351  and 
formatted  into  the  linear  option.   At  this  point  the  user 


11 


selects  the  (R)eceive  option  on  the  Z248,  presses  the  "PRINT 
DATA"  prompt  on  the  screen  of  the  351  and  enters  a  filename 
on  the  Z248. 

The  final  version  of  the  data  transfer  program  was 
compiled  using  Microsoft's  BASIC  compiler  Version  l.O.   By 
eliminating  the  writting  of  the  data  to  the  screen  in 
addition  to  the  file  and  using  the  compiled  program  vice  the 
BASIC  interpreter  a  file  which  will  use  60K  of  disk  space 
can  now  be  transferred  in  under  5  minutes  as  opposed  to  the 
20  minutes  it  took  with  SMARTCOM  II.   A  listing  of  the 
program  is  included  in  Appendix  B. 


12 


III.   CORROS I ON^RAIES 

A.   THEORETICAL  EQUATIONS 

In  general  the  corrosion  behavior  of  a  metal  can  be 
described  by  the  use  of  the  Nernst  equation,  activation 
potential ,  concentration  potential  and  mixed  potential 
theory. 

The  Nernst  equation  is  a  modification  of  the  Gibb's 
Free  Energy  equation  and  expresses  the  tendency  of  a 
corrosion  reaction  to  proceed  in  terms  of  EMF  based  on  the 
activities  of  the  products  and  reactants  of  the  1/2  cell. 


E   =  E'='  -  RTlog.(a) 
nF 


E   =  Electrode  oxidiation  potential 

(0  for  reduction) 
E'^  =  Equilibrium  oxidiation  potential 

(0  for  reduction) 
R   =  Gas  constant 
T   =  Temperature  '^'K 

n   =  Number  of  electrons  taking  part 
F   =  The  Farday 
Q   =  Activities  of  products/activities  of 

reactants 

The  potential  of  the  cell  is  the  algebraic  sum  of 
oxidation  and  reduction  1/2  cells.  ..  i, 
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Activation  polarization  describes  the  current  density 
which  exists  when  the  cell  is  polarized. 

i*r-ie3€dxc=  =  i«,exp( +ZFaE/RT) 
ic:«t:r^o.jj.<=  =  icexp(-ZF(l-a)E/RT) 


io  =  equilibrium  exchange  current  density 

Z  =  Valance  of  the  active  species 

E  =  Polariztion  about  the  equilibrium  potential 

a  =  fraction  of  E  which  applies  to  i«n«aciA<= 

R,T,F  as  above 


By  manipulating  the  above  equation  the  activation 
polarization  -*«  can  be  written  as: 

■7,  =  0  logao  (irt/io) 
0   =  2.3RT/aZF 

Although  activation  polarization  occurs  for  both  the 
anodic  and  cathodic  reactions  the  cathodic  reaction  can  also 
be  effected  by  concentration  polarization.  When  this  occurs 
the  reduction  reaction  occurs  at  such  a  rate  as  to  decrease 
the  concentration  of  ions  in  the  solution  near  the  surface 
of  the  metal.   This  results  in  a  change  in  the  potential  of 
the  1/2  cell  reaction  since  the  effective  activity  in  the 
Nernst  equation  changes. 
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The  concentration  polariztion   =  is  expressed  as: 


'*?  c  =  2 ■  3RT  loQio  (1  -id-^ii.) 
nF 


id  =  operating  current  density 
x,_    =  limiting  current  density 


The  limiting  current  density  is  inversely  proportional 
to  the  thickness  of  the  diffusion  layer.   This  thickness  is 
dependent  on  various  geometric  and  environmental  factors  and 
can  often  only  be  determined  experimentally.   The  effect  of 
concentration  polarization  can  be  minimized  by  agitating  the 
solution,  thus  decreasing  the  layer  thickness,  or,  as  in  the 
case  of  oxygen  reduction,  by  heating  the  solution  to  reduce 
the  concentration  of  the  dissolved  gas.   Both  of  these  means 
of  dealing  with  concentration  polarization  may  significantly 
alter  the  desired  corrosion  environment  and  distort  the 
correlation  between  the  experimental  results  and  the  actual 
in— service  corrosion  rates. 

In  the  absence  of  concentration  polarization  a  dynamic 
polarization  diagram  of  a  cell  reaction  would  appear  as  in 
Figure  1.   The  slope  of  the  anodic  and  cathodic  curves 
d0/d(logxo  id)  would  be  given  by  B  as  described  in  the 
activation  polarization  expression. 

Most  cells  are  affected  by  concentration  polarization 
such  that  the  cathodic  polarization  '^  =  iyj±    ~7*    ^T*=*  ^^^ 
effect  of  this  is  shown  in  Figure  2.   It  is  obvious  from 
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this  that  loqio  (id)  approaches  Logxo  (ii_)  the  linear 
behavior  of  the  polarization  curve  is  distorted  and  in  cases 
where  ii_  and  i^^  are  not  sufficiently  separate,  a  linear 
region  may  not  exist. 

Mixed  potential  theory  concerns  concerns  the  immersion 
of  a  single  piece  of  metal  into  an  electrolyte  and  the 
subsequent  reactions.   The  theory  consists  of  two  simple 
hypotheses. 

1.  Any  electrochemical  reaction  can  be  divic'ed  into  two  or 
more  partial  oxidiation  and  reduction  reactions. 

2.  There  can  be  no  net  accumulation  of  electrical  charge 
during  an  electrochemical  reaction.  [Ref.  12  :  pp  314] 

It  is  these  hypotheses  that  allows  for  the 
determination  of  two  important  parameters  necessary  for 
estimating  corrosion  rates.   When  the  cathodic  branch  of  the 
1/2  cell  with  the  highest  equilibrium  potential  intersects 
the  anodic  branch  of  the  other  1/2  cell  the  corrosion 
potential  0c=c»r-r-  and  the  corrosion  current  density  iccsi-i-  can 
be  used  to  estimate  corrosion  rates. 

When  the  corrosion  system  consists  of  only  a  single 
metal  in  contact  with  a  solution  which  contains  a  single 
reduction/oxidation  reaction  the  overall  reaction  is  the 
linear  sum  of  the  reduction  currents  and  the  oxidation 
currents. 
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If  the  solution  contains  more  than  one 
reduction/oxidation  reaction  the  overall  reaction  may 
proceed  as  the  linear  combination  provided  that  the 
equilibrium  exchange  current  densities  of  the  secondary 
reduction  reactions  are  within  one  order  of  magnitude  of  the 
main  reduction  reaction.   The  secondary  oxidiation  reactions 
proceed  in  the  same  manner.   As  a  result  the  secondary  redox 
reactions  affect  either,  both,  or  none  of  the  branches  of 
the  overall  reduction/oxidiation  reaction.  CRef.  13:  pp. 
199] 

If  concentration  polarization  is  not  present  both  the 
anodic  and  cathodic  polarization  curves  Mould  consist  of  a 
of  horizontial  lines,  representing  the  current  added  by  the 
reaction  at  that  potential,  and  slanted  lines  of  some  siiope 
0  representing  the  activation  polarization  of  all  the 
oxidation  or  reduction  reactions  which  are    influential  at 
that  potential.   In  the  presence  of  concentration 
polarization  the  possibility  of  introducing  numerous  non- 
linearities  into  the  cathodic  branch  is  distinct.   As  the 
overall  reduction  reaction  is  the  sum  of  all  the  influential 
reduction  currents  any  region  could  consist  of  zones  where 
the  activation  polarization  of  the  ith  reaction  dominates 
the  concentration  polarization  of  the  (i— l)th  reaction. 
This  can  make  the  potentiodynamic  curve  extremely  difficult 
to  analyze. 
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B.   DETERMINATION  OF  CORROSION  RATES 

As  previously  described  the  presence  of  numerous 
electrodes  and/or  metals  in  the  corrosion  system  can 
significantly  alter  the  shape  of  the  polarization  curve. 
Despite  this  if  it  was  possible  to  design  a  current 
measurement  system  that  was  capable  of  measuring  only  the 
anodic  and  cathodic  currents  the  determination  of  iccar-i- 
could  still  be  simply  made  since  log^o  i«=c»t-f-  would  be  at  the 
intersection  of  the  anodic  and  cathodic  branches  of  the 
polarization  diagram. 

In  the  real  world  polarization  curves  take  on  a  'Y' 
shape  as  a  result  measuring  instrument  limitations  and  the 
complexity  of  the  corrosion  environment.   0c:oi^i-  can  be 
measured  directly  since  it  is  the  equilibrium  potential  of 
the  system.   As  the  potential  is  scanned  towards  0cc»-r-  the 
net  current  (i.rtociic  -  ic=««^r^c3c«j.c= )  is  measured.   Only  when 
the  potential  is  displaced  far  enough  from  0cc»r-f-  does  one  of 
the  two  currents  overwhelm  the  other  and  does  typical  Tafel 
behavior  take  place. 

The  presence  of  numerous  1/2  cell  reactions, 
concentration  polarization  and  the  need  for  one  of  the 
currents  to  overwhelm  the  other  can  all  combine  to  make  the 
determination  of  ice».-r-  difficult. 

If  concentration  polarization  is  not  present  at  some  60 
from  0COI-I-  the  polarization  diagram  should  adopt  true  linear 
behavior.   The  amount  of  60  at  which  the  curve  should  adopt 
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linear  behavior  is  in  itself  open  to  disscussion.   Ailor 
states  that  80  should  be  at  least  50mV  from  0c:or-r-  [Ref.  14: 
pp.  199],   Once  the  linear  regions  are    found  two  methods  can 
be  employed  to  determine  i^a.^^-       The  simplest  is  to  draw  two 
straight  lines  from  the  linear  regions  towards  0cor-.-.   The 
two  lines  should  intersect  at  ic=or-f   A  second  method  is  to 
measure  the  slopes  0  (Tafel  constants)  of  the  linear  regions 
of  the  anodic  and  cathodic  branches.  These  values  can  be 
used  in  conjunction  with  data  determined  from  a  polarization 
resistance  test  to  determine  the  value  of  icor-r-  by  use  of 
the  equation: 

icor-i-  — _0«=ft« 

2.3(B«+Bc)PR 

f 
If  concentration  polarization  is  predominant  over  the 

cathodic  region  of  the  polarization  curve,  the  so— called 

'knee'  method  is  often  employed.  In  this  method  a  tangent  is 

drawn  on  the  knee  of  the  cathodic  branch  such  that  the 

angles  formed  on  both  sides  of  the  tangent  with  the  cathodic 

branch  are    equal.  When  this  has  been  accomplished  either  of 

the  methods  discussed  in  the  previous  paragraph  may  be  used. 

Many  other  methods  have  been  developed  to  determine 

iccar-r-  or  to  determine  corrosion  rates  directly  from 

polarization  data.  Princeton  Applied  Research  uses  a  Chi- 

squared  minimization  technique.  This  technique  is  included 

in  the  software  package  which  operates  their  Model  351 

corrosion  measurement  system  as  PARCalc.  By  by  their  own 
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admission,  since  the  program  is  terminated  based  on  a  change 
in  the  average  value  of  the  square  of  the  relative 
deviations  between  two  successive  iterations,  it  is 
necessary  to  determine  the  Tafel  region  before  starting  the 
process.  [Ref.  15:  pp  IX-IO] 

LeRoy  [Ref.  16: pp  1006-1012D  postulated  the 
determination  of  Tafel  Slopes  based  on  polarization 
resistance  techniques.  Without  repeating  the  entire 
discussion  it  is  sufficient  to  say  his  method  of  determining 
Tafel  slopes  based  on  expected  values  from  activation  and 
concentration  polarization  expressions,  based  on  the 
assumption  that  the  parameters  necessary  to  evaluate  those 
expressions  were  available,  did  not  accurately  describe  the 
observed  values.  [Ref.  17:  pp  1988-1989] 

It  would  also  seem  possible  to  solve  the  overvoltage 
expression  for  the  terms  in  question.  Recalling  the 
overvoltage  expression: 

0  -  0c=e,r-.-  =  Bc=logj.o(ic/ic)  +  2.3RTloqxo(l  -  id/ii_) 

nF 

The  problem  is  that  0^,  ie»  and  i^  all  represent  unknowns 

and  the  ii_  term  is  located  in  such  a  manner  as  to  complicate 

the  solution. 
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The  loQio  (1  -  irt/ii_)  term  can  be  replaced  by  the  series 
expansion : 

logxo(l+X)  =  2.3C  X  -  X =  +  X^  -  X^  ...  +(-l)"*^X"l 

2     is     4  n 

where  X  =  — ic*/ii_ 
The  final  equation  to  be  solved  would  be: 

0  -  0cr«=,.-r-  =  0c:logxo(ic^)-Pc=lagxo(ioi)+2.3C5:(-i)"-^X"] 

n 

The  difficulty  with  this  approach  is  that  a 

sufficient  number  of  terms  must  be  carried  in  the  series  in 

order  to  accurately  represent  the  logxo  term.   This  causes 

two  severe  problems.   The  first  is  n+2  data  points  must  be 

used  to  solve  the  system  of  equations.   When  one  considers 

that  id  is  normally  on  the  order  of  uA/cm^,  exponentiating 

it  to  the  nth  degree  can  easily  result  in  a  matrix  to  be 

solved  involving  coefficients  of  50  or  60  orders  of 

magnitude.   The  end  result  is  that  while  the  system  can  be 

solved  such  that  the  coefficients  determined  will  result  in 

the  accurate  representation  of  the  actual  overvoltage,  the 

coefficients  themselves  have  been  influenced  by  round— off 

error  and  catastropic  cancellation,  and  represent  merely 

numbers. 

The  other  problem  exists  with  attempting  a  solution  of 

this  form.   As  discussed  previously  with  the  possibility  of 

mixed  electrodes,  the  associated  redox  processes  may  or  may 

not  effect  the  overall  behavior  of  the  system.   Even  if  the 
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limiting  current  densities  ii_  were  all  known,  it  would  still 
be  necessary  to  determine  their  effect  on  the  overall 
process  in  order  to  include  them  in  the  system  of  equations 
to  be  solved.   As  a  result,  for  m  different  electrodes,  m 
experiments  involving  just  the  electrode  and  the  metal  would 
have  to  be  carried  out  just  to  determine  their  influence. 
Consider  the  use  of  such  a  method  with  synthetic  seawater. 
At  least  10  different  chemical  substances  could  effect  the 
overall  reaction  in  addition  to  the  major  metal  oxidiation 
and  hydrogen  reduction  reactions.   It  should  be  obvious  that 
if  the  limiting  current  densities  are  not  known  the  size  of 
the  system  of  equations  to  be  solved  could  easily  exceed  the 
number  of  available  data  points.   More  importantly,  one 
would  have  to  assess  the  advantages  of  performing  a 
polarization  in  a  mixed  medium  if  the  same  test  has  to  be 
performed  separately  on  the  individual  components  first. 

C.   CALCULATION  OF  TAFEL  CONSTANTS 

In  order  to  determine  the  location,  if  one  existed,  of 
the  anodic  and  cathodic  Tafel  regions  two  numerical 
dif f erentation  techniques  were  employed.  These  two  methods 
were  the  Four  Point  Central  Difference  method  and  the  use  of 
a  Cubic  Spline  interpolating  polynomial.   A  graphical 
display  of  the  derivatives  allowed  the  user  to  determine  the 
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region  of  linearity.   Once  this  region  was  determined  a 
first  order  linear  regression  is  performed  on  the  original 
data  within  the  region  to  determine  the  Tafel  slopes. 
1 •   Four  Point  Central  Difference  Method 

When  the  data  points  are    equally  spaced  the  Central 
difference  method  may  be  used  to  approximate  the  derivatives 
at  a  point.   The  use  of  this  method  can  result  in 
considerable  time  savings  over  the  cubic  spline  method  as 
the  derivative  can  be  calculated  without  having  to  resort  to 
the  solving  of  a  system  of  equations  for  each  point  in 
question.   The  method  is  not  as  accurate  as  the  cubic  spline 
method  and  graphical  representations  of  the  derivatives  will 
often  contain  more  "noise"  than  those  derived  from  the  cubic 
spline  method. 

The  four  point  Central  Difference  formula  for 
computing  derivatives  is: 

d0 =  _12h 

ddogxo  ic.)     Q{  fx    -     f-x)    +    ( /-2  -  Vs) 

Where: 

h   =  0A  —  0I-H 

fr,    =  logio(ici)  at  the  Ith  location  from 
logj.o(  id)  I 

The  Central  Difference  method  has  an  error  term 

associated  with  it.   In  this  case  the  error  term  is  on  the 

order  of  l/h'*.   Although  this  may  appear  to  be  a  substantial 

term,  it  should  be  remembered  it  is  a  passible  error,  the 
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method  is  only  being  used  to  determine  the  Tafel  region  and 
when  used  the  Tafel  constant  is  inverted  so  the  error  term 
would  actually  be  h'* . 

2 .   Cubic  Spline Interpolating  Polynomia.l 

Wfien  trying  to  numerically  determine  derivatives  of 
experimental  data  which  may  or  may  not  be  egually  spaced  the 
use  of  interpolating  polynomials  is  often  suggusted.  One  of 
the  most  accurate  methods  is  to  pass  a  cubic  spline  thru  the 
data  points  and  then  to  differentiate  the  resulting  eguation 
to  arrive  at  the  derivative  at  a  desired  print.  CRef.  18: 
pp.  242] 

The  advantage  of  using  a  cubic  spline  is  that  in 
deriving  the  system  of  eguations  to  arrive  at  the  constants 
at,b,Cf    and  d   which  form  the  eguation: 

y  =  ai(x  -  Xi)=  +  ^i(x  -  Xi)=  +  c, (x  -  Xj)  +  dx 

the  system  is  formed  in  such  a  manner  that  the  value 
of  the  function,  first  derivative  and  second  derivative  ^re 
the  same  for  the  pair  of  cubics  which  join  at  each  point- 

This  reguires  that: 

y'  =  3ai(x  -  Xi)=  +  2^1  (x  -  Xi)  +  Cx 
and       y'  '  =  6ax  ( x  -  Xj)  +  2bx 

In  order  to  simplify  the  the  mathematics  involved 
the  eguations  for  the  system  are    written  in  terms  of  the 
second  derivatives  of  the  interpolating  cubic  Si. 
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The  resulting  interpolating  cubic  will  then  take  on 
the  form  of : 

hi_xSi-i+2Si (hi_i+Sx )+hiSi^i=6(yi-i-yi )-  (yi-yi-x) 

hi  hi  — J. 

where  hi  =  x  -  Xi 

Solving  the  system  leads  to  the  following  operations 

to  determine  the  values  of  ^x f bx , Cx    and  dx • 

ax  =  (Si^o.  -  Si)/6hi 

bx    =  Si/2 

cx    =  ((yx-x  -  yi)/hx)  -  ((2hxSx  +  hxSx*x)/6) 

di.    =    Yx 

The  resulting  system  of  n  -  2  equations  in  Sx 
involves  n  pairs  of  data  points  in  order  to  generate  the 
required  number  of  equations.   To  arrive  at  the  two 
additional  equations  to  solve  for  Sx  and  Sr, ,  constraints  Are 
specified  which  pertain  to  the  conditions  at  the  ends  of  the 
curves.   The  three  choices  for  the  end  conditions  arez 

1.  Sj.  =  Sr,  =  O,  which  implies  that  the  end  cubics 
approach  linearity  at  their  extremities- 

2.  Sx  =  Sz,  Sr.  =  Sp,_x.   This  assumes  that  the  end  cubics 
approach  parabolas  at  their  extremities. 

3.  Sx  is  a  linear  extrapolation  of  Sx  and  S^.  S„  is  a 
linear  extrapolation  of  S„—z    and  Sn— x - 
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In  this  case  since  the  interest  is  in  determining 
the  derivatives  d0/d(logio  ia)     pointwise,  and  not  in 
developing  an  interpolating  polynomial,  it  was  decided  to 
use  the  end  conditions  of  Sj.  =  S„  =  O,  and  to  use  enough 
data  points  as  to  force  the  inconsistences  caused  by  the 
assumed  end  conditions  away  from  the  point  in  question. 
Four   data  points  on  each  side  of  the  point  in  question  were 
used  in  determining  the  the  derivative  at  the  point  in 
question.   This  resulted  in  the  solving  of  a  7  x  7  system  of 
equations  for  each  data  point.   By  locating  the  point  in 
question  in  the  center  of  the  data  set,  the  accuracy  of  the 
estimate  involved  was  improved,  and  the  mathematics  involved 
in  determining  the  derivative  were  simply  solving  for  Cs.. 
4 .   Linear  Regression 

Once  the  anodic  and  cathodic  Tafel  regions  have  been 
determined,  the  data  points  which  comprise  these  regions  Are 
used  in  a  least  squares  curve  fitting  routine  to  determine 
an  equation  of  the  form: 

0  =  ^(logxo  (i.*))  +  b, 

where  0  represents  the  Tafel  slope.  The  values  of  a 
and  b  are    found  by  salving  the  two  following  equations 
simul taneous I y . 

^(logxo(ic*)i  )=  +  /3S(  logio(ici)i  )  =  Sdog^oC  id  )  i  )  (0x  ) 

^(  logxo(id)i)  ^    tN  =  Z(0x) 

where  N  is  the  number  of  paints  in  the  Tafel  region. 


26 


5 .   Considerations  Affect  ing  _.Me  t  hod Se  l_ec  t  i  on 

When  the  PAR  351  Corrosion  Measurement  system  scans 
the  potential  during  a  Potentiodynamic  test  the  current  is 
measured  at  intervals  of  2mV  or  0,5mV  of  potential, 
depending  on  the  relative  change  in  the  current.   The 
sampling  rate  is  normally  0.5mV  when  the  current  shows  a 
substantial  change  and  2mV  when  it  does  not.   This  usually 
results  in  a  sampling  rate  of  0.5mV  in  the  Tafel  regions, 
although  severe  concentration  polarization  in  the  vicinity 
of  i«=or-r-  may  result  in  a  limited  number  of  points  having 
0.5mV  intervals  on  the  cathodic  branch.   Two  other  aspects 
may  also  effect  the  potential  intervals  at  which  the  current 
appears  to  have  been  sampled. 

Mhen  data  is  sent  out  the  serial  port,  the  351 
normally  uses  a  format  similar  to  the  FORTRAN  Ga.b  type. 
Not  counting  leading  zeros,  four  significant  figures  are 
used.  This  results  in  a  problem  in  the  recorded  potential 
values.   When  the  potential  is  greater  than  or  equal  to  1.0 
volts,  or  less  than  or  equal  to  -1.0  volts,  only  three 
decimal  places  Are    carried.   As  a  result,  when  the  values  of 
0e:or-t-  fall  in  this  range  and  the  scan  rate  is  0.5mV,  the 
recorded  data  seems  to  reflect  that  the  potential  has  not 
changed  for  two  successive  current  measurements. 

Since  the  potential  is  being  scanned  it  is  rather 
obvious  when  this  occurs  that  the  recorded  data  does  not 
accurately  represent  the  experimental  data.   Since  this 
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pattern  is  easily  recognized,  the  array  containing  the  data 
can  be  manipulated  tc  more  accurately  reflect  the 
experimental  data. 

Occasionally  the  potentials  recorded  display  no 
regular  pattern  and  the  data,  although  guestionable  in 
nature,  must  be  used  as  recorded. 

When  calculating  the  derivatives,  the  anodic  and 
cathodic  branches  are    performed  separately.   The 
calculations  begin  at  the  third  data  point  away  from  Ecoi-r- 
when  the  Central  Difference  method  is  used,  and  at  the  fifth 
when  using  the  cubic  spline  method.   The  calculations 
continue  as  long  as  the  data  points  continue  to  suggest  a 
slope  of  of  the  expected  sign  (+/-).   This  automatically 
excludes  the  passivating  region  of  the  anodic  section,  as  xt 
should  since  it  is  not  part  of  the  Tafel  region,  and 
excludes  the  possibility  of  the  program  terminating 
execution  if  it  were  to  calculate  an  infinite  slope.   In 
general  this  strategy  worked  as  expected  and  provided  enough 
data  points  to  properly  examine  the  Tafel  regions.   The  only 
experiment  in  which  it  failed  was  for  the  SONOSTON  sample 
tested  in  3.5yC  NaCl  solution.   As  can  be  seen  by  examining 
the  Potentiodynamic  curve  of  this  sample  (Figure  3) ,  the 
cathodic  branch  exhibited  behavior  throughout  which  caused 
either  method  to  terminate  execution  without  generating  any 
derivatives.   The  reason  for  this  behavior  is  unknown. 
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I V .   GENERAL  OVERVIEW  OF  THE  PROGRAM 

The  program  was  developed  so  as  to  allOM  for  its  use 
without  the  supposition  that  the  user  was  familar  with  the 
FORTRAN  langauge  or  with  personal  computers  in  general. 
Numerous  efforts  were  made  to  minimze  the  possibilty  of  the 
program  prematurely  terminating  execution  in  the  event  of 
improper  input  or  through  floating-point  operations 
resulting  in  invalid  expressions.   By  takxng  advantage  of 
the  color  capability  of  the  monitor  screen,  colors  are 
changed  with  respect  to  the  type  of  input  expected,  and  to 
display  messages  resulting  from  either  input  errors,  data 
file  format  errors  or  a  failure  to  properly  follow  the  steps 
necessary  for  proper  program  execution. 

The  program  was  written  using  Microsoft  FORTRAN 
Optimizing  Compiler  version  4.O.   Graphics  were  achieved 
using  Microcompatibles,  Inc.'s  GRAFMATIC  library.   Plotter 
support  was  achieved  using  Microcompatibles,  Inc.'s 
PLOTMATIC  library. 

A  summary  of  the  main  program  and  its  subroutines  is 
given  below.   The  main  program  and  its  called  subroutines 
are  listed  in  Appendix  C.   The  subroutines  are  presented  in 
the  order  in  which  they  appear  in  the  main  program. 
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A.  ST ART. FOR 

This  is  the  main  program.   The  name  of  the  first  data 
file  to  be  used  is  entered.   Following  calls  to  DATAIN, 
CHECK  and  DAT AM  a  menu  consisting  of  6  options  is  presented: 

1.  Plot  the  potentiodynamic  curve 

2.  Generate  Tafel  Slopes  using  cubic  spline 

3.  Generate  Tafel  slopes  using  central  difference 

4.  Use  another  data  file 

5.  Overlay  two  potiodynamic  curves 

6.  Exit 

Option  1  generates  calls  to  GRAPHl  and  PLOT 

Option  2  calls  SLOPE  1 .   Following  this  a  first  orde*- 
linear  regression  of  the  anodic  and  cathodic  branches  using 
locations  within  the  original  array  determined  during  a  call 
to  GRAPHl  is  performed  and  the  results  along  with  original 
potentiodynamic  curve  ^re    passed  to  GRAPHl. 

Option  3  calls  SLOPE.   Upon  return  from  SLOPE  this 
option  proceeds  as  Option  2 

Option  4  returns  the  user  to  the  3rd  executable  line  of 
START . 

Option  5  calls  DATAIN,  CHECK  and  DATAM  while  retaining 
the  data  from  the  first  file. 

Option  6  obviously  exits  the  program. 

B.  DATA IN. FOR 

This  opens  the  input  file  and,  using  a  formatted  read 
statement,  generates  the  array  containing  potential  and 
current  density.   Additionally  this  subroutine  counts  the 
total  number  of  lines  in  the  file.   This  is  done  to  ensure 
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that  formatting  changes  within  the  data  file  which  were  not 
corrected  during  the  data  file  transfer  will  not  have  not 
caused  the  array  generation  to  prematurely  terminate.   If  a 
formatting  error  has  been  detected,  an  error  message  is 
displayed  in  red  and  the  line  number  within  the  file  is 
identified.   This  allows  the  user  to  exit  the  program  and 
correct  the  error  before  continuing. 

C.  DAT AM. FOR 

This  converts  the  absolute  valua  of  the  current  density 
to  its  logxo  value  and  determines  the  position  of  Ecor-i-  with 
the  array. 

D.  CHECK. FOR 

This  subroutine  determines  the  number  of  data  points  in 
the  anodic  and  cathodic  branches  to  be  used  in  the 
determining  the  derivatives  when  using  the  central 
difference  method.   This  is  done  by  checking  log  current 
density  for  values  which  might  result  in  the  calculation  of 
an  infinite  slope  (normally  the  onset  of  pitting  or 
passivation)  and  for  equal  increments  in  the  recorded 
potentials.   Occasionally  the  recorded  potentials  will  not 
reflect  the  polarization  which  has  taken  place.   This  occurs 
when  the  absolute  value  of  the  recorded  potential  is  equal 
to  or  greater  than  1.0.   In  these  ranges  the  O.SmV  voltage 
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increment  falls  outside  of  the  total  width  of  the  field  and 
is  not  recorded.  As  this  is  a  recognizable  pattern,  at  the 
users  option  the  problem  can  be  corrected. 

E .  SLOPE . FOR 

This  is  the  first  subroutine  called  when  using  the 
Central  Difference  method.   It  assigns  graph  titles,  axes 
labels  and  calls  DATADEL.FOR. 

F.  DATADEL-FOR 

Using  the  points  determined  in  CHECK,  this  calculates 
the  numerical  derivatives  using  the  Central  Difference 
method . 

G.  SLOPE 1. FOR 

This  assigns  graph  titles  and  labels.   It  also 
determines  the  number  of  points,  MI,  to  be  used  on  both 
branches  in  the  cubic  spline  method  by  checking  only  for 
values  of  log^o  current  density  which  could  result  in  the 
calculation  of  slopes  of  infinite  value.   It  calls  CSPLIN. 
separately  for  the  anodic  and  cathodic  branches.   MI  anodic 
is  not  necessarily  egual  to  MI  cathodic. 

H.   CSPLIN. FOR 

For  MI  number  of  points,  this  generates  the  7x7  matrix 
and  the   7x1  vector  of  values  described  in  the  cubic  spline 
section  MI-9  times.   As  each  matrix  is  formed  it  is  solved 
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by  a  call  to  LINSYl.FOR.   Upon  return  from  LINSYl  the  7x1 
vector  contains  the  solution  which  is  used  to  calculate  the 
derivative. 

I.   LINSYl.FOR 

Using  partial  pivoting  this  subroutine  solves  the 
augmented  7x8  array    by  Gaussian  Elimination.   The  original 
7x1  vector  contains  the  solution  CRef.  19], 

J.   GRAPHl.FOR 

This  is  the  graphics  output  of  the  selected  main  menu 
option.   Through  the  use  of  the  numeric  keypad  as  softkeys 
the  user  can  scale  the  axes  as  desired.   The  option  to 
return  to  the  original  graph  or  the  previous  display  always 
exists.   When  plotting  the  results  of  either  numeric 
dif f erentation  method,  only  those  points  which  consecutively 
fall  within  the  desired  region  are    plotted.   Once  the  axes 
have  been  scaled  such  that  graph  contains  only  the  linear 
region  of  the  anodic  or  cathodic  branches  the  position  of 
the  start  and  stop  points  within  the  original  array        ' 
containing  the  derivatives  ar&    saved.   Since,  when 
calculating  the  derivatives,  the  data  point  containing  Ecorr 
is  always  used  as  a  reference,  ie.  (depending  on  the  method) 
the  first  derivative  for  each  initial  branch  is  either  3  or 
5  points  away  from  Ece»r-r-,  the  start  and  stop  points  also 
represent  the  position  of  the  Tafel  regions  within  the 
original  potential  and  loglO  current  density  array.   When 
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the  graph  is  not  of  the  numerical  derivatives,  any  data 
point  which  falls  within  the  desired  axes  is  plotted.   If 
desired  GRAPH!  calls  PLOT, FOR  for  a  hard  copy  output  of  what 
is  currently  displayed  on  the  screen - 

K.   PLOT. FOR 

This  plots  the  graph  displayed  on  the  screen.   Although 
many  of  the  arguments  passed  from  GRAPHl  are  used  directly 
as  passed,  this  is  not  a  'screen  dump'.   The  resolution  of 
the  plct".  .'r  is  much  higher  than  that  of  the  screen.   As  a 
result  what  appears  to  be  a  one  pixel  jump  on  the  screen 
will  not  plot  as  such.   Two  sets  of  tabulated  data  may  be 
plotted  on  each  set  of  axes.   At  the  users  discretion 
different  color  pens  may  be  used  to  plot  the  axes,  first 
output,  second  output,  and  labels. 

L.   CONFIGURATION 

In  its  final  form,  the  program  has  been  compilied  and 
linked  into  an  executable  file  which  runs  on  a  Zenith  248 
computer  which  is  equipped  with  a  80286  microprocessor,  an 
80287  Numeric  Processor  Extension,  a  EGA  graphics  card  with 
64K  video  RAM,  and  a  color  monitor  which  supports  the  EGA 
color  mode.   The  plotter  used  is  an  Houston  Instruments  DMP— 
40  series  digital  plotter. 

The  plotter  library  is  plotter  series  specific.  As  a 
result  of  this,  the  use  of  a  plotter  other  than  a  Houston 
Instrument's  DMP— 29  or  greater  would  require  obtaining 
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another  library  from  flicrocompatibles  Inc.   As  many  of  the 
subroutines  contained  within  PLOTMATICS  Are    similar  to  those 
called  within  GRAFMATICS,  it  can  probably  be  safely  assumed 
that  little  if  any  editing  of  the  subroutine  PLOT. FOR  would 
be  required. 

To  use  the  program  with  a  video  configuration  other  than 
that  previously  described  would  require  editing  of  the 
subroutines  PLOT. FOR,  GRAPHl.FOR,  DATAIN.FOR,  START. FOR  and 
CHECK. FOR.   DATAIN,  START  and  CHECK  would  require  editing 
only  if  the  video  configuration  was  monochrome.   PLOT  and 
GRAPHl  would  require  editing  for  any  graphics  configuration 
other  than  EGA. 
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V.   RESULTS  AND  DISCUSSION 

A.   SYNTHETIC  SEAWATER 

Table  1  below  compares  the  Tafel  constants  determined  by 
using  the  new  program  with  those  obtained  by  Escue  CRef.  93 • 

TABLE  1 

TAFEL  CONSTANTS  OBTAINED  FROM 
SYNTHETIC  SEAWATER  EXPOSURE 


MATERIAL 

0. 

0-» 

0c 

0=* 

Ti-50Ni 

0.1400 

0.174 

0.1155 

0.132 

1020  Steel 

0.0758 

0.068 

0.1481 

0.114 

304  Steel 

0.0241 

0.178 

0.0207 

7075  Al 

0.3790 

0.087 

0.1584 

0.062 

Fe-Cr-Mo 

0.1775 

0.121 

0.1226 

0.141 

Fe-Cr-Al 

0.1B40 

0.180 

0.1369 

0.151 

Cu-Mn-Al-Fe-Ni 

0.0399 

0.048 

0.1579 

0.120 

Cu-Mn-Al 

0 . 0633 

0.035 

0.0652 

0.045 

Cu-Zn-Al 

0.0804 

0.062 

0.0785 

0.075 

630  Bronze 

O . 0403 

0.031 

0.0724 

0.041 

»  From  Escue  [Ref .  9] 

Table  2  below  compares  the  current  densities  and 
related  corrosion  rates  of  the  samples  tested  in  synthetic 
seawater.   The  first  line  for  each  material  represents  the 
current  density  as  determined  by  the  intersection  of  the 
Tafel  lines.   In  many  of  the  cases  the  anodic  and  cathodic 
Tafel  lines  did  not  intersect  at  a  single  point  and  are 
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joined  by  a  horizontal  line.   In  these  cases  icor-r-  was 

determined  by  taking  the  average  of  the  anodic  and  cathodic 

ice-,-  at  the  point  where  the  lines  intersected  a  horizontal 

line  form  by  0co,-i-- 

The  second  line  represents  the  results  obtained  by 

combining  the  Tafel  slopes  with  the  polarization  resistance 

data. 

TABLE  2 

CURRENT  DENSITY  AND  CORROSION  RATES 


MATERIAL 

•Icor-r- 

■Lfzcyr-r- 

CR 

CR* 

uA/cm^ 

tiA/cm= 

mpy 

mpy 

Ti  -  50y.Ni 

1.245 

1.734 

0.573 

0.799 

4.182 

4.960 

1.822 

2.160 

1020  Steel 

4.212 

3.837 

1.981 

1.805 

20.325 

17.201 

9.580 

8.110 

7075  Al 

1.413 

0.421 

0.591 

0.176 

3.067 

0.994 

1.281 

0.415 

304  Steel 

0.106 

0.728 

0.048 

0.329 

0.157 

1.190 

0.071 

0 .  537 

Fe-Cr-Mo 

0.750 

o.ao4 

0.354 

0.380 

2.294 

2.060 

1.085 

0.974 

Fe-Cr-Al 

0.822 

0.830 

0.362 

0.365 

1.183 

1.240 

0.520 

0.545 

Cu-Mn-Al-Fe- 

Ni   8.241 

8.854 

4.165 

4.070 

9.662 

10.400 

4.882 

5.260 

Cu-Mn-Al 

2.089 

1.119 

1.111 

0.595 

1.459 

0.894 

0.784 

0.481c 

Cu-Zn-Al 

2.239 

1.897 

1.136 

0.963 

3.780 

3.230 

1.919 

1.640 

630  Bronze 

1.762 

1.442 

0.901 

0.737 

6.878 

4.690 

3.550 

2.420 

*  From  Escue 

CRef.  9] 
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1  •    I  i  -  5P/CN  i  (JIN  I LG ) 

These  results  compared  favorably  with  those  obtained 
by  Escue.   Figure  4   shows  the  derivatives  as  calculated  by 
the  cubic  spline  method.   As  can  be  seen  in  this  figure  no 
true  linear  regions  Are    distinguishable  in  either  of  the  two 
branches.   When  this  occurs  the  "knee"  method  can  be 
employed  by  examining  the  branches  for  a  range  where  the 
derivatives  oscillate  about  a  linear  midpoint.   In  this  case 
both  the  anodic  and  cathodic  branches  showed  this 
oscillation  in  the  log  current  density  range  of  -5.28  to  - 
5.51.   Figure  5  shows  the  resultant  Tafel  lines  plotted  on 
the  original  potentiodynamic  curve.   Although  the  curves 
intersect  yJccar-r-  at  what  could  be  considered  to  be  a  single 
point  the  icc-r-  determined  in  this  manner  differs  greatly 
from  that  calculated  by  using  the  Tafel  constants  and  the 
polarization  resistance  data. 

2 .   1020  Carbon  Steel  (51020LG) 

Figures  6  and  7  show  the  results  as  obtained  for  the 
1020  carbon  steel  sample.   Of  all  of  the  samples  tested  this 
one  exhibited  a  cathodic  branch  with  the  least  tendency  to 
behave  in  a  linear  fashion.   If  there  was  any  linear 
behavior  in  the  cathodic  branch  it  occured  at  a  log^o  ico.-.- 
of  about  -5.2.   The  anodic  branch  displayed  a  much  longer 
lasting  region  around  log,.o  icc-r-    =    -    4.3.   This  sample  also 
showed  the  greatest  disparity  in  the  value  of  icot-t- 
dependinq  on  the  manner  in  which  it  was  calculated. 
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3 .  7075  -  Alum  i  nun*  ( AL  707  5Lj3 J. 

In  this  case  the  central  difference  method  showed  a 
linear  region  in  both  the  anodic  and  cathodic  branches.   As 
can  be  seen  in  Figure  8  the  linear  regions  exist  at  the 
pronounced  peaks  in  the  derivatives-   In  this  case  the  Tatel 
lines  did  intersect  at  a  single  point  as  seen  in  Figure  9. 

4 .  304  Stainless  Steel  (5ST304LG) 

The  calculated  corrosion  rates  and  the  respective 
Tafel  constants  differed  greatly  from  those  obtained  by 
Escue.   Although  the  calculated  corrosion  rates  correlated 
better  to  the  actual  weight  loss  methods,  the  need  to 
terminate  the  calculating  of  derivatives  before  the  onset  of 
passivity  may  have  caused  a  linear  region  in  the  anodic 
branch  to  be  omitted.   These  points  can  be  seen  by  the  sharp 
peak  on  the  anodic  branch  on  Figure  lO,  and  the 
corresponding  point  on  the  paten tiodynamic  curve  Figure  11. 
The  Tafel  lines  intersected  at  a  single  value  of  log^o  icor-i- 
but  this  time  there  was  a  good  correlation  between  the 
calculated  and  graphical  iccsi-i-- 
5  -   Fe-Cr-Mo  (VACROILG) 

Again  the  presence  of  any  truly  linear  regions  is 
questionable,  although  both  branches.  Figure  12,  show  a 
range  where  the  deviations  from  linear  behavior  ^re    minimal. 
Figure  13  once  again  shows  a  theoretical  intersection  of  the 
Tafel  lines. 
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^-   Fe-Cr-Al  (_VACRq2LG) 

This  alloy  showed  behavior  similar  to  the  Fe-Cr-Mo 
alloy  described  above.   What  may  be  considered  a  truly 
linear  region  in  the  cathodic  branch  is  shown  in  Fiqure  14 
just  to  the  left  of  the  sharp  peak.   The  overlay  of  the 
Tafel  lines  on  the  potentiodynamic  curve  is  shown  in  Fiqure 
15.   The  Tafel  constants  derived  show  a  closer  relationship 
to  those  calculated  by  Escue,  than  did  a  comparision  of  the 
Fe— Cr — Mo  sample. 

7  .   Mn-Al -Fe-Ni  _(SgNO_S JLG 

SONOSTON,  Figures  16  and  17,  displayed  a  linear 
anodic  branch  and,  as  did  most  of  the  alloys,  a  cathodic 
branch  effected  by  concentration  polarization. 

8  .   Cu-Mn-Al  _  (.INCRLG ) 

The  result  here  again  showed  a  linear  anodic  branch 
but  in  this  case  the  concentration  polarization  did  not 
overwhelm  the  cathodic  branch  in  such  a  manner  as  to  distort 
the  region  where  activation  polarization  is  predominant. 
Figure  18  shows  the  graph  of  all  the  derivatives  and  Figures 
19  and  20  show  the  graphs  of  the  derivatives  in  what  was 
decided  to  be  the  Tafel  region.   Figure  21  shows  the 
original  potentiodynamic  and  the  Tafel  lines. 

9  .   Cu-  Z  n -Al_.iPLALC:LGl 

Experinced  similar  behavior  as  the  Cu— Mn— Al  sample. 
Once  again  a  linear  region  was  very  apparent  in  the  anodic 
branch.  Figures  22  and  23. 
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10.  630  Series  Bronze  (BRNZLG) 

This  was  the  only  sample  of  those  tested  in  the 
synthetic  seawater  where  the  Tafel  regions  for  the  anodic 
and  cathodic  branches  were  found  in  the  same  region  of  log 
iccar-r--   As  Seen  in  Figure  24  concentration  polarization 
again  guickly  dominates  the  cathodic  branch.  Figure  25 
displays  the  overlay. 
B.   3.57.  NaCl  SOLUTION 

In  a  similar  fashion  Table  3  lists  the  Tafel  constants 

and  Table  4  the  resultant  current  densities  obtained  from 

the  samples  tested  in  the  3.57C  NaCl  solution.   In  general 

the  current  results  show  a  stronger  correlation  with  the 

corrosion  rates  obtained  by  Akthar  [Ref.  lO]  than  did  the 

comparision  of  the  current  results  with  those  obtained  by 

Escue  CRef.  9]  for  the  synthetic  seawater. 

TABLE  3 

TAFEL  CONSTANTS  OBTAINED  FROM 
3.5  -/.  NaCl  SOLUTION  EXPOSURE 


MATERIAL 

0. 

B-* 

0= 

Bo* 

304  Steel 

0.1259 

0.4117 

0.0951 

0.117 

7075  Al 

0.0116 

O.OIO 

0.0383 

2.36 

Fe-Cr-Mo 

0.2821 

0.234 

0.1053 

0.110 

Fe-Cr-Al 

0.1203 

0.214 

0.1292 

0.1285 

Cu-Mn-Al 

O . 0082 

0.0137 

O . 0337 

1.379 

630  Bronze 

0.0524 

0.570 

0.1472 

0.321 

»  From  Akthar  [Ref.  10] 
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TABLE  4 
CURRENT  DENSITY  AND  CORROSION  RATES 


MATERIAL 


7075  Al 


304  Steel 


Fe-Cr-Mo 


Fe-Cr-Al 


Cu-Mn-Al 


630  Bronze 


1  ccai-r- 

-•■Gca  r-r- 

CR 

UA/Clfl^ 

uA/cm= 

mpy 

6.960 

7.406 

0.048 

10.460 

11.700 

1.842 

0.052 

O.IOO 

0.026 

0.012 

0.020 

0.006 

0.043 

0.040 

0.019 

0.609 

0.613 

0.274 

0.493 

0.660 

0.234 

0.916 

1.180 

0.428 

2.864 

6.062 

1.454 

5.251 

lo.aoo 

2.666 

0.093 

0.129 

0.048 

3.570 

4.480 

1.842 

CR* 
mpy 

2.965 
4.685 

0.050 
0.010 

0.018 
0.276 

0.313 
0,551 

3.078 
5.483 

0.066 
2.312 


♦From  Akthar  [Ref.  10] 


1 -  7075  Aluminum  (AL7075) 

In  this  case  the  effect  of  concentration 
polarization  made  any  estimate  of  the  cathodic  Tafel 
constant  extremely  difficult.   As  seen  in  Figures  26  and  27 
concentration  polarization  completely  overwhelmed  the 
cathodic  branch.   As  a  result  only  the  "knee"  method  cuold 
be  used,  but  its  accuarcy  is  questionable.   The  anodic 
branch  displayed  typical  linear  behavior. 

2 -  304  Stainless  Steel  (SST304) 

The  cathodic  branch  of  this  sample  did  display  a 
linear  region  located  between  the  sharp  peaks  in  Figure  28. 
The  anodic  branch  displayed  numerous  attempts  to  passivify 
as  seen  in  Figure  29. 
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3 .   Fe-Cr-Mo  (VCMOPD) 

Similar  behavior  to  the  304  Stainless  steel  above. 
The  bandwidth  exhibited  by  the  cathodic  branch  near  loga.o 
icot-r-  of  about  -  6.5  Figure  30,  typifies  the  knee  in  the 
curve.   Again  numerous  attempts  to  passivify  limited  the 
range  of  data  which  could  be  use  in  the  anodic  branch  Figure 
31. 

^'       Fe-Cr-Al  (VACAL2) 

Both  the  anodic  and  cathodic  regions.  Figure  32, 
displayed  a  region  where  the  behavior  could  be  considered  to 
be  linear.   Although  the  Fe— CR— Mo  and  Fe— Cr — Al  samples 
tested  in  synthetic  seawater  displayed  similar  behavior,  the 
Fe-Cr — Al  sample  tested  in  the  3.5%  NaCl  solution  did  not 
exhibit  the  tendency  to  passivify.  Figure  33,  as  did  the  Fe- 
Cr — Mo  sample. 

5.   Cu-Mn-Al  (INCRMTE) 

This  sample  best  typifies  the  problem  associated 
with  concentration  polarization.   On  Figure  34  a  linear 
region  can  be  seen  in  the  cathodic  branch  at  a  logxo  icof-i- 
of  about  -5.3.   This  region  is  followed  by  a  rapid  decrease 
in  the  slope  of  the  curve  indicative  of  concentration 
polarization.   When  the  original  potent iodynamic  curve  is 
viewed  thru  out  its  entire  polarization  range.  Figure  35, 
very  little  can  be  determined  due  to  the  great  changes  in 
current  density  on  the  anodic  branch.   When  the  range  in  the 
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cathodic  branch  displaying  the  linear  bevaior  is  scaled 
properly  a  textbook  example  of  concentration  polarization  is 
seen.  Figure  36.   The  overlay  is  shown  in  Figure  37. 
6 .   630  Series  Bronze  (BRZ630PD) 

As  seen  on  Figures  38  and  39  this  sample  experienced 
similar  behavior  to  the  one  above.   All  though  the  linear 
region  directly  before  concentration  polarization  is  not  as 
well  defined  the  cathodic  behavior  is  dominated  first  by 
activation  polarization  and  then  by  concentration 
polar i  zation . 

C.   COMPARISION  OF  CURRENT  DENSITIES 

In  many  of  the  samples  tested  the  corrosion  rates  as 
determined  by  the  intersection  of  the  Tafel  lines  were 
significantly  lower  than  those  obtained  by  the  use  of  the 
Tafel  Slopes  and  linear  polarization  data.   As  previously 
stated  the  Tafel  regions  should  begin  to  dominant  the 
polarization  curve  at  an  polarization  increment  of  about 
50mV  from  Ecoi-r-.   This  implies  that  a  linear  region  should 
not  be  present  at  increments  less  than  this.   With  these 
tests  a  linear  region  was  often  found  within  this  zone. 
Because  of  the  general  shape  of  a  potentiodynamic 
polarization  curve  the  resultafit  lines  drawn  from  this 
region  will  always  have  a  lower  slope  and  thus  intersect  at 
a  lower  value  of  i.   The  disparity  underscores  the  need  to 
use  the  average  of  icor-r-   as  determined  by  the  intersection 
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of  the  Tafel  lines  and  that  calculated  by  the  Tafel  slopes 
and  linear  polarization  data  when  estimating  corrosion 
rates . 

Another  interesting  feature  was  noticed  when  conducting 
some  data  manipulation  of  the  potentiodynamic  polarizations 
analyzed.   When  the  absolute  value  of  (0  -  0cor-r-)  was 
plotted  against  log^o  iej  in  many  of  the  cases  the  anodic  and 
cathodic  branches  began  to  deviate  significantly  at  a  log  id 
corresponding  to  the  intersection  of  the  two  Tafel  lines. 
This  may  be  coinicidental  but  also  may  warrant  further 
discussion.   When  a  linear  polarization  resistance  test  is 
conducted  the  potential  is  scanned  with  a  range  of  0^^^^-    +/- 
25  mV.   The  polarization  resistance  is  determined  by 
calculating  the  slope  of  the  linear  region  about  i  =  O 
A/cm2.   If  a  linear  polarization  curve  is  plotted  in  the 
same  manner  as  above,  while  neglecting  the  singularity  at  i 
=  O  A/cm^,  the  deviation  discussed  would  determine  the  end 
of  the  linear  region.   Returning  to  the  equation  used  to 
determine  icon-  using  the  Tafel  constants  and  the 
polarization  resistance  it  can  be  shown  that  in  order  to  the 
polarization  resistance  to  be  a  constant  that  0«0c  = 
C*(0«+0c:),  where  C  is  a  constant  of  proportionality.   The 
current  being  measured  during  a  potentiodynamic  polarization 
is  equal  to  the  difference  between  the  anodic  and  cathodic 
currents.   Only  when  one  of  the  currents  begins  to  overwhelm 
the  other  would  there  be  a  substantial  change  in  the 
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measured  current  with  respect  to  a  given  overvoltage.   This 
Mould  correlate  to  the  deviation  between  the  anodic  and 
cathodic  branches  on  a  plot  of  the  absolute  value  of  (0  - 
0c:c3r-r-)  vs  1  og  ( irf )  .   It  a  1  SO  should  correspond  to  the   end 
of  the  linear  region  on  the  linear  polarization  curve  since 
at  this  point  the  constant  of  proportionality  would  no 
longer  be  a  constant  as  one  of  the  two  branches  is  beginning 
to  overwhelm  the  other.   It  is  postulated  that  a  technique 
of  this  method  could  be  used  to  determine  the  minumium 
overvoltage  at  which  the  search  for  a  Tafel  constant  could 
begin.   It  is  recognized  that  this  does  not  account  for  the 
possibility  of  completely  symetric  anodic  and  cathodic 
branches  but  if  this  was  the  case  there  would  be  no  problem 
in  resorting  to  the  +/-  50mV  rule  proposed  by  Ailor. 
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V I .   CONCLUS I ON  AND  RECDMMENDAT I ONS 

The  FORTRAN  based  graphics  program  proved  to  be 
extremely  valuable  in  determining  the  existence  of  linear 
Tafel  regions  particularly  in  those  cases  where 
concentration  polarization  dominated  the  cathodic  branch. 

With  the  ability  to  transfer  the  experimental  data  to  a 
programable  computer  the  opportunity  exists  to  further 
analyze  the  data.   This  could  allow  for  the  determination  of 
which  electrodes  Are    effecting  the  main  redox  reaction  when 
the  environment  consists  of  a  solution  of  mixed  electrodes. 
In  order  to  accomplish  this  it  is  recommended  that  a  single 
high  damping  alloy  be  subject  to  individual  poten tiodynamic 
polarizations   containing  a  single  component,  at  its  normal 
concentration,  of  synthetic  seawater.   Through  the  use  of 
graphical  technigues  it  should  be  possible  to  isolate  the 
influential  redox  reactions. 

In  order  to  accomplish  this  it  will  be  necessary  to 
conduct  all  the  experiments  in  solutions  which  contain  the 
same  amount,  or  at  least  a  monitored  amount,  of  dissolved 
oxygen.   This  is  one  factor  which  was  not  considered  in  the 
previous  two  experiments  and  may  have  made  the  analysis  more 
difficult  than  necessary. 
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In  lieu  of  this  approach  an  effort  should  be  made  to 
determine  the  limiting  current  densities  causing  the  severe 
concentration  polarization. 
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APPENDIX    B 
PAR    351    TO    2-248    DATA    TRANSFER    PROGRAM 


10  J  =  O 

20  SCREEN  O 

30  aX.OR  4,1,1 

40  KEY  CFF:CLS:CLDSE 

50  DEFINT  A-Z 

60  LOCATE  25,1 

70  PRINT  STT?INB*(60,"    ") 

80  FALSE  =  0:TFLJE=NDT  FALSE 

90  XCFF*=<>«*(19):X0N!fc=a-R*(17) 

lOO  LOCATE  25,20:FRINT   "PAR  351   TO  Z-24«  DATA   TRANSFER  PROGRAM" 


110  LOCATE  3,20:FRI^^" 
120  LOCATE  5,20:PRINT 
130  LOCATE   7,20:PRINT 
140  LOCATE  a,20:FRINT 
150  LOCATE  9,20:F«INr 
160  LOCATE   10,20:PRI^fT 
170  LOCATE   ll,20:PRIIVfT 
180  LOCATE   12,20:PRINT 
190  LOCAIE   13,20:PRINT 
200  Ca_DR   14,1,1 
210  LOCATE   15,20:RRINT 


"VOU  MUST  PROPEH_Y  CO^FIGLR    TVE  PAR  351." 
"AFTER  SETTirsG  TVE    Tir-E  AM)  DATE  GN  TFE  351" 
"FfESS     SYSTEM  MPNAGER' " 
"PRESS     COrvFIQJR  SYSTEM" 
"SET    ' PROTOCOL'    TO  PRINT" 

"SET      PARITY'    TO  ND^E" 

"SET    'BPUD  W^TE      TO  2400" 

•SET     STCP  BIT'    TO  1" 

'SET    'ICRD  LEMJm'    TO  8" 


"PFESS  ANY  KEY  TO  CQNTIIN4JE" 
GOTO  220 


220  GO*  =   IJVJ<EY*:IF  GD*= 

230  CLS 

240  CGLDR  4,1,1 

250  LOCATE  25,20:PRIIvrr   "PAR  351   TO   Z-248  DATA  TRANSFER  PROGRAM" 

260  LOCATE  3,20:PRINr   "PRESS     MAIN  rtMJ'    CN  TFE  351" 

270  LOCATE  4,20:PRINT   "PFESS    'RECAL1_/DISF1_AY  EXFtRIfENT' " 

280  LOCATE  5, 20:  PR  I  r^  "PRESS     DISPLAY  EXPERirtNT '  CN  TVE     COPY  FROM  DUK' 

SECTION" 

"SELECT  TVE  EXPERirCNT  TO  BE  TfJANSFEP^ED" 

"PRESS     PLDT/FCFJIAT'" 

"EN5LFE  TVE   X  AXIS   IS   IN  T>t  LINEAR  SCALE" 


'290  LOCATE  6,20:PRINT 

300  LOCATE  7,20:PRINT 

310  LOCATE  a,20:PRIMT 

320  COLGR  2,1,1 

330  LOCATE  9,20:PRINfr 

340  LOCATE   10,20:RRIMT 

FORMAT    " 

350  COLOR  14,1,1 

360  LOCATE   14,20:PRINT 

370  GO*=IN<EY$:IF  GD*= 

380  CLS 

390  COLOR  4,1,1 

400  LOCATE  25,20:F«INT 

410  g^EED*  =   "2400" 

420  COTFIL*  ="CDri2:"+SreED»-^",N,8,l,CS,DS" 

430  OPEN  COTFIL*  PG  #1 


"IF  NOT    IN  A  LINEAR  DISPLAY  TVEN  PRESS" 
"FDFJiAT  DISFl^Y,    LINEAR- ,'FtVIEW  DISPLAY' ,    FIXJT 


FTESS  ANY  KEY  TO  CONTIMJE" 
"   GOTO  370 


PAR  351   TO  Z-248  DATA  TRANSFER  PROGRAM* 
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AAO  PRINT   ttl,    XCFF*; 

450  iJH=N   "SOPNz"   RjR  OLfTRJT  AS  #2 

460  UXATE   10,25:PRINT   "(R)€3ceive  a   file  OR    (E)xit   " 

470  LiXAIH   L1,37:IM=LIT   TXRX* 

480   IF      (TXRX*<>"R")    fifO    (TXRX1i<>"E")    THEN  460 

490   IF  TXRXt="E"   TVEN  1350 

50O  CUB 

510  aUDR  4,1,1 

520  LOCATE  25,20:RRirn"   "PfiR  351   TD  Z-248  DftTA  TfWNS^R  R«3GRAM" 

530  COLOR  15,1,1 

540  LOCATE   13,25:FT?I^^"   "FTESS     PRINT  DATA'    ON  T>t  351   TVEN" 

550  LOCATE   15,30:PRINT   "PRESS  PNY  KEY  TO  COTiTIMJE" 

560  GO*  =    IM<EY*:IF  GD*=""   GOTO  560 

570  CIS 

580  COLOR  4,1,1 

590  LOCATE  25,20:PRINr   "PAR  351   TO  Z-248  DATA  TRW^EFER  FR0GFVY1" 

600  LOCATE  10,20:PRINT   "THIS  PRDERAM  USES  A  Wfl  DISK  TO  S3PEED 

EXECUTION. " 

610  COLOR  2,1,1 

620  LOCATE  13,20:PRINT   "DO  NOT  ENTER   TT€  Dia<  DRIVE  OR  T>C  FIUE  TYPE" 

630  LOCATE   14,20:PRINT   "BUT  OM.Y  TT-E  Nf^Tt  OF  THE  OLm^LTT  FILE" 

640  COLOR   14,1,1 

650  LOCATE   16,20:PRIMT   "ENTER  Tl-E  OLITPUT  FILE  N«t" 

660  LOCATE   18,25:  INPLTT  FIL* 

670  DSK*="D:":DeKFIL*=DSK$^FIL* 

6B0  OPEti   "D-.m^"   FOR  OUTPUT  AS  #3 

690  N  =  O 

700  CLS 

710  PRINT  «1,    XONS; 

720   IF  L0C(l)=O  TFEN  GDSLB  820 

730   IF  LDC(1)>128  TVGM  PAUBE=TRUE:  PR  I  NT  #1,X0FF* 

740  A*=INRjr*(l,»l) 

750   IF    (A*=a-R»(10))   OR    (A*=a-F»(13) )    TFEN  B«=a-R»(13)  :J=J+1 

760   IF    {A*<>a^«i(13))    AM)   (A»<>CHR*(10) )    T>€N  PRINT  #3, A*; 

770   IF    (J   >   3)   AND   (B*=CHFWi(13) )    AM)    (A*=CW«i(13) )    T1-€N  PRINT 

#3,a^«i(13) 

780   IF    (J    >   3)   AND   (B»=C>fWi(13) )    AM)   (A*=a-Rt(13) )    Tl-EN  B*="    ":N=N«-1 

790   IF  L0C(1)>0  Tl-EN  730 

800   IF  PAUSE  TVEN  PAUGE=FALSE:PRINT   »1,X0N*; 

810  GOTO   720 

820  FOR   1=1   TD  5000 

830   IF  LUC(1)<>0  TTEN   I   =  9999 

840  ^EXT    I 

850    IF    I>9999  TFEN  RETURN 

860  CUDGE  #3:CLS 

870  F*="D:"+f» 

880  OPEN  "D:T1^"  FOR  IISPUT  PG   #4 

890  OPEN  D6KFIL*  FOR  OUTFUT  fG   #3 

900  FOR  I  =  1  TO  N 

910  X*  =  I^F1JT*(44,#4) 

920  Y»  =  IM^UT*(2,#4) 
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930  B*  =  MICHi(X*,ia,9) 

940  CU  =  MID1i(B*,i,l) 

9S0  D*  =  Mim(B*,2,l) 

960  E*  =  MID*(B*,3,1) 

970  F*  =  MID*(B*,4,1) 

9B0  G*  =  MID*(B*,5,1) 

990  l-Hi  =  MID*(3*,6,1) 

lOOO   m  =  mD1i(B*,7,l) 

1010  J*  =  MID*(B*,8,1) 

1020  KU  =  MID*(B*,9,1) 

103O  FX*  =  "-0,000000" 

1040   IF    (E*  =   "E")    TVCN  MID*(X*,18,9)    =  FX* 

(F*  =  "E")    TVEN  MID*(X*,18,9)    =  FX* 

(6*="E")    IT-EN  MID*(X*,18,9)    =  FX* 

(H*="E")    TT-EN  liirWi(X1i,18,9)    =  FX* 

(H*="E")    TFEN  MID*(X*,18,9)    =  FX* 

(I*="E")    TFEN  MID*(X*,18,9)=FX* 

(J*="E")    THEN  MID*(X*,18,9)=FX* 

(K*="E")    Tl-EN  «ID*(X*,ia,9)=FX* 
1120  F2*  =   "E-12   " 
1130  A*  =  MID*(X*,42,1) 
11^  A1*=«ID*(X*,43,1) 
llbO  A2*=tiID*(X*,34,l) 
1160  A3*  =  A*+A1* 

1170   IF    (A3*  =   "12")   AM)   (A2*  =   "-" )    TFEN  f1ID*(X*,39,5)    =  F2* 
1100  P*  =  MID*(X*,1,1) 
1190   IF    (P*  =   "P")    TFEN  PRINT#3,X* 
1200  ^EXT    I 
1210  CLOSE  #3:CLS 
1220  J   =  O 
1230  COLCR  1,5,4 
1240  BEEP 

IZiX)  LOCATE  12,30:PRINT   "ttttlTlf^tGFER  Cai=LETE*«««" 
1260  LOCATE  13,30:PRINT   "FTCSS  ANY  KEY  TO  CGNTIfsLE" 
1270  G0*=1M<EY*:IF  BD*=""  GOTO  1270 
12BO  PRllsrr  #1,    XOFF* 
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IF 

1060 

IF 
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IF 

1090 

IF 
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12^?0 

CLUdt 

#1 

130O 

CLOSE 

#2 

1310 

CLOSE 

#3 

1320 

n.lTF 

#4 

1330 

m  rF 

580 

,1 

1340 

buiu  : 

1350 

END 

92 


APPENDIX    C 
CORROSION    PROGRAM 


THIS    IS   1>E  MftIN  PR0GRW1.      WITHIN   IT    Tl-E  USER  ENTERS  T>t  FILE 
NPrt  WHICH  ft/a-E  Wishes  TO  PNflLYZE.      FOLLDWI^e  THIS  A  rtNU 
APPEARS  ALLOWING  TVE  USER  TO  SELECT  Tl-E  DE5IFED  OPTION. 


TTE  FOLLOWING  VV^IBLES  ARE  DEFI^ED 

CDATA    :    TTE  2  DirCNSICNPL  AFRAY  FDR  TVE.  DATA  FILE  C0NTA^NI^6 

POTENTIAL  AND  CURRENT   DE>eiTY 
CDATAl:   AS  ABO^  BLTT  OM_Y  LEED  W-EN  C^ER_AYI^e  TWO  CURVES 
MDATA   :   A  TWO  Dir-ENSIONAL  ARFV^Y  OF  CDATA  BUT  NOW   IN  TEPWS  OF 
POTENTIAL  AM)  LQGIO  CURRENT  DENSITY 
A  OTE  DirtNSIOMAL  ARRAY  CC^^^AINI^G  TFE   X  VALUES 
OF   THE  FIRST  Of^VE  TO  HE  GRflH^D  OR  FUJI  ILD 
A  OrvE  DirCNSIONPt.  ARFViY  CONTAINING  TTE  Y  VALUES 


XDAT 

YDAT 

XDATl 

YDATl 

XI 

X2 

XY 

Yl 

NS 

B2 

Bl 

M 

I 

N 

01 

Nl 


J 

Jl 

02 

J2 

SI 


OF  TFE  FIRST  CURSE  TO  EE  GRAPFED  OR  F1.0TTED 

A  C^E  DII-EN5IGNPL  ARFiAY  C0NTAINIM5  THE   X  VALUES 

OF  TFE  SECDM)  CURVE  TO  BE  GRAPHED  OR  F1_0TTED 

A  OhE  DirENBION«_  ARRAY  CONTAINING  TFE  Y  VALUES 

OF   TFE  SECDM)  CLRVE  TO  BE  GRAFTED  OR  PLOTTED 

ThE  SLM  0=^  TFE  LOG  DJTENT  DE^eITIES.    USED   IN  TFE 

LI^EAR  REGRESSION  FDR  DETEFJ1INING  TFE  TAFEL  CD^eTANTS 

TFE  SUM  OF  TFE  SQUARES  OF  TFE  LOG  CURRENT  DE^eITY 

INITIALLY  USED   IN  TFE  SPf-E  tlftf^R  A  XI    BUT  ONCE  IhE. 

TAFEL  CONSTANT  HAS  BEEN  DETEFWirCD  USED  AS  TFE  LOG 

CURFENT  DENSITY  VWLUE   IN  CREATIM3  TFE  TAFEL  LI^ES 

ThE  SLM  OF  TFE  PRODUCTS  OF  LOG  CUF»esfT  DENSITY  AM) 

POTENTIAL.   USED   IN  TFE  SATt  MPM>ER  AS  XI   AM)  X2 

TFE  SLfl  OF  TFE  POTENFIAL.    USED  AS  XI,X2,XY 

T>E  NUhCER  OF  POINTS  USED  TO  6E^ERATE  TFE  TAFEL  SLa=E 

TFE  TAFEL  SLOPE  AG  CPLOULATED   IN  TFE  REGRESSION 

TFE   INTEFCEFT  OF  TFE  TAFEL  LIISE 

rtfMJ  RESPOTEE 

DO  LOOP  COUNTER 

TFE  M-MBER  OF  POINTS   IN   TFE  DATA  FILE  BEING  FCPO 

TFE  NUhOER  OF  POINFS   IN  TFE  FIRST  DATA  FILE,    UNCHPNGIM5 

TFE  NUrCER  OF  POINTS   IN  TFE  FIRST  DATA  FILE.    UEED    IN 

TFE  SUBROUTI^E  SLOPEI   AM)  CAN  BE  RETLRsED  WITH  A 

DIFFERENT  VflLUE 

TFE  LOCATION  OF  ECOF«  WITHIN  fCATA 

SATE  AS  J   BUT  PASSED  AS  AN  ARGUE^ENT  TO  SLDPEl 

SATt  AS  01   BUT  USED  OM_Y  WFEN  0^EF8_AYING  DATA  FILES 

SATE  AG  J   BUT  USED  DM_Y  WFEN  D^ERLAYING  DATA  FILES 

TFE  POSITION  WITHIN  rt)ATA  WFERE  TFE  REaUIRE^ENTS  OF 

TFE  CENTRAL  DIFFERNCE  fETHOD  FAIL.   CATHDDIC  BRANCH 
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«  S2  :    SAft  P6  SI   BLfT  FOR   THE  flMDDIC  BRANCH 

»  S3  :    TVE   TOTAL  I^JffiER  OF  POINTS   IN  TVE  CATHDDIC  BRAfCH 

«  iOHICH  tEEJ  TVE  FEGUIRB-ENTS  FDR  LEINB  TVE  CENTFJAL 

DIFFEFENCE  r-EHHOD 
S4  :    SATE  AS  S3  BUT  FOR  TFE  AMJDIC  BRANCH 

Al  :    TVE  POSITION  WITHIN   XDATl   AM}  YDATl  WHICH  MAR<S  TVE 

:    BEGIM^ING  OF  TVE  TAFEL  REGION 
A2  :    SAME  fiG  AL   BJT  TVE  Ehf)  OF   TFE  FEBION 

A3  :    SERJEB  TTE  SAft  PURPOSE  AS  Al  BUT  WITH  RESPECT  TO 

riDATA 
A4  :   SATE  AS  A3  BUT  AS  A2 

CI  :    SA«-E  AS  Al  BUT  WITHIN  XDAT  ATO  YDAT   (CATVODIC) 

C2  :    SATE  AS  CI   BUT  TVE  EM)  CF   TVE  REGION 

C3  :    SERVES  TVE  SA^E  PURPOSE  PG  CI   BUT  WITH  RESPECT  TO 

rf>ATA 
C4  :    SAft  AS  C3  BUT  AS  C2 

NB  :    SET  TD  O.   PASSED  TO  GRAFHl  AM)  PLOT    ID  AVOID  TVE  USE 

OF   TFE  EhPTY   XDATl  A^^  YDATl  ARRAYS  WFEN  TIE  ARE  NOT 

USED 
NPME      :    TVE  NATE  OF  TVE  MAIN  DATA  FILE   IN  USE 
NAI-E2   :    TVE  NATE  OF  TFE  SECOND  FILE  W-EN  0«.€R_AYING 
TITLE2:    TVE  SECDM)  VE/^INB  FOR  TVE  GRAPH  AM)  FLOT.    CHPM3ES 

DEPEMENT  ON  TVE  OPTION  FF\  H'Abl) 
TITIJE3:    RESERVED  FOR  FUTURE  USE    IF  DE5IFCD 
Dl  :    CHARACTER  STRING  USED  TD  FCFW  TnLE2 

02  :    TFE  CHARCACTER  REFf^SENTATION  OF  B2  WHEN  B2  REPRESENT 

TVE  CATVCDIC  TAFEL  SLOPE 
D3  :   CHARACTER  STRING  USED  TO  F0FJ1  TITLE2 

D4  :    SATE  AS  D2  BUT  FOR  TVE  ANODIC  TAFEL  SLOPE 

LIST      :   CHARACTER  STRIM3  » P»^  FCR  MESSAGES  AM)  FROTPTS 


C 
C 
C 
C 
C 
C 
C 
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REAL*4  XDAT (2000), YDAT (2000), XDATl (2000), 
+YDAT1(2000) 

REAL*8  M)ATA(2000,2),rt)ATAl(2000,2),Xl,X2,NS,XY,XI,Yl 
+  ,CDArA(200O,2)  ,CDATA1( 2000,2) 

INTEGER»2  I , J, N,A1,A2, CI, C2, 01, NI ,J1, SI, 32,53,34,02, J2,NB 

INTEGER»2  A3,A4,C3,C4 

a4ARACTER»8  NP#t,NArE2,    TITLE2«32,YNP^E«5,LIST*50,TITLE3«50 

CHARACTER»8  D2,D4,D1«5,D3«6 

TITLE2  =  '  • 

YNAft  =  '  ' 


TVE  SCREEN  IS  QLEAFED  AM)  THE  FILE  NATE  IS  ENTERED.  TVE  USER 
CAN  ALSO  EXIT  TVE  PROtfJAM  OR   GET  A  LIST  OF  EXISTING  DATA  FILES. 
THIS  PROr-PT  APPEARS  IN  FU^R_E  WHICH  IS  TVE  COLOR  USED  Ut«^EVER 
A  (3HARACTER  RESPONSE  IS  EXPECTED  TD  A  PRDTFT.  IF  A  FILE  NATE  IS 
ENTEFED  WHICH  DOES  NDT  EXIST  AN  ERROR  l-ESSAGE  WILL  APPEPR  IN  RED. 

CALL  QSM)DE(3) 
CALL  QCL£AR(0,5) 
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15  LIS!    =    '     ' 

CPLL  a:MDV(15,12) 

U«ITE(6,805)   LIST 

Cf^LL  GO-DV(25,i2) 

LIST  =     ENTER  7VE  FILE  Nf^tE.  TO  BE  LISED' 

WRITE(6,a05)    LIST 

CPLL  QCrt]V(20,ll) 

LIST  =     A  BLft4<:   TO  EXIT,   OR   "DIR"   FOR  A  DIRECTGRY" 

WRnE(6,805)   LIST 

CPU-  QCJt:W(36,10) 

READ(6,800)    Nfift 

CALL  QD-CV(40,9) 

IF    iNPt^    .EQ.    '  '  )    T>-€N 

GOTO  9CX) 
E^O^F 
IF    (NflTE    -EQ.      DIR)    TTEN 

CALL  Gn_EAR(0,5) 

CALL  QCrf3/(l,23) 

PAUSE     ENTER   "DIR  «.    /W"   OR  A  BLAM<  TO  CCNTINUE' 

CALL  QLJSCFL(23,6,1,1,79,0,5) 

GOTO  15 
ENDIF 

CALL  DATAIN(CDATA,N,NPrt,IER) 
IF    (lER    .EQ.    1)    TVEN 

GOTO  lO 
EJOIF 

CALL  DATAM(CDATA,M)ATA,N,J) 
CALL  a-ECK(rDATA,N,J,Sl,S2,S3,S4) 

20  CALL  aSMDDE(3) 
CALL  QCLEAR(0,2) 
Jl  =  J 

N  =  N-1 

01  =  N 

Nl  =  N 
C 

C  TFE  MAIN  rtNU   IS  DISR_AYED   IN  GREEN.      Tl-E  COLOR  GREEN   IS  USED 

C  U*E>eV^R  AN   INTEGER   IS  EXPECTED  AS  THE  FESPCN5E  TO  A  PRCH=T. 

C  IF  TVE    IIM=UT   IS  MJT  AN   INTEGER  OR  A  CHflRftCTER  AN  tF¥SJR  MESSAGE 

C  IS  DISH_AYED   IN  RED. 

C 

21  CALL  Qa-CV(15,22) 
LIST  =  '  • 
WRITE(6,805)  LIST 
CALL  QCMDV(15,ia) 

LIST  =   1.  PLOT  ORIGINAL  POTENT lODYNPMIC  CLf^E' 

WRITE(6,805)  LIST 

CALL  QD1DV(15,17) 

LIST  =      '2.    GEI^eViTE  TAFEL  SLOPES  USING  CUBIC  SFLIfvE  rtTTOD' 

WRnE(6,805)    LIST 

CALL  GCMCW(15,16) 

LIST  =      '3.   GE^EFV^^E  TAFEL  SLCFES  USING  CENTRflL  DIFFEREMZE  fETTDD' 
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UFITE(6,a05)    LIST 
C^^-L  aCM3/(15,i5) 

LIST  =      '4,    LEE  flNOn-ER  DATA  FILE' 
URI7E(6,a05)    LIST 

LIST  =  '5.  OVERLAY  TWO  POTENT  I  ODYNfiMIC  CURVES' 
CALL  QDt]V(15,14) 
URI7E(6,a05)  LIST 
CPLL.  QCriCW(15,13) 
LIST  =  '6.  EXIT' 
WRITE(6,805)  LIST 
CALL  aCM^'(29,12) 
LIST  =  'ENTER  1,2,3,4,5  OR  6' 
WRITE(6,805)  LIST 
CfiLL  aCMDV(39,10) 

REPO(6,810,BRR  =  22,IQGTAT  =  J4)    M 
C«_L  QCrOV(39,9) 
22  IF    ((J4    .^E.   O)    .OR.    (M   .ETT.   6))    TVEN 

CALL  QEMDDE(3) 
CALL  QaLEPR(0,4) 
LIST  =    '     ' 
CALL  GCMDV(20,i7) 
WRITE(6,805)   LIST 

LIST  =    'YOU  DID  NOT  ENTER  AN   INTEGER  OF' 
CALL  QCM]^/(20,16) 
WRnE(6,805)   LIST 
LIST  =    '\/flLLJE  1,2,3,4  OR  5.' 
CALL  QCI1DV(20,15) 
WRITE(6,805)   LIST 

LIST  =     PFESS  ANY  KEY  TO  CONTINUE' 
CPLL  QCMDV(20,13) 
WRITE(6,805)   LIST 
CflLL  aCMDV(39,12) 
CflLL  QIM<EY(I3,I4) 
GOTO  20 
E^DIF 
C 

C  OPTION  1   HAS  BEEN  SELECTED.   THIS  6E^ERATES  A  GFWPH  AM)  OPTIONAL 

C  PLOT  CF  TKt  ORIGINflL  POTENTIODYNfiMIC  CURVE. 

C 

IF    (M   .EQ.    1)    TVEN 
DO  25   I   =   1,01 

XDAT(I)    =  M)ATA(I,2) 
YDAT(I)    =  M)ATA(I,1) 
25  CCNTIMJE 

Nl   =  N 
^B  =  O 

YNfirt  =    'VO-TS' 
TITLE2  =    '     ' 
TITLE3='     ' 

CALL  GRflPHl(XDAT,YDAT,01,Jl,Nflft,XDATl,YDATl,re,YNfift,TITLE2 
-»-  ,A1,A2,C1,C2,TI7LE3) 
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G07D  20 
EM)IF 

C 

C  OPTICN  2  H«S  BEEN  SELECTED.    THIS  6E^E^^TES  TVC  EERIVATIVES  O- 

C  THE  PQTENTIODYNPMIC  CURVE  US1^G  TVE  CUBIC  SR_I^E  ftTrHDD. 

C 

IF  (M  .EQ.  2)  THEN 
Al  =  O 
A2  =  O 
CI  =  O 
C2  =  O 
30  Jl  =  J 
Nl  =  N 

CPLL  SLCFEl(M)ATA,Nl,Jl,NPtt,Al,A2,Cl,C2) 
C 

C  STILL  WITHIN  CFTICW  2.    TVE  PF03RAM  D-ECKS  TO  ENSLFE  THAT  TVE 

C  T«^EL  FEBICN5  FCR  BOTTH  THE  ANODIC  ATsO  CATFCDIC  REGIONS  HAVE 

C  BEEN  SELECTED. 

C 

IF    ((CI    .GE.    C2)    .OR.    (Al    .GE.    A2) )    TFEN 
CALL  QSrDDE(3) 
CALL  QCLEAR(0,4) 
LIST  =    '     ' 
CALL  GOCV(15,20) 
WRITE(6,805)   LIST 

LIST  =     YOU  DID  NOT  SELECT  BOTH  AN  ANODIC 
CALL  aCMDV(15,i9) 
U«ITE(6,805)   LIST 

LIST  =     AND  CATHDDIC  EPPfCH  OR  ELSE  TVE  POINTS' 
CALL  GCM:V(15,18) 
WRITE(6,805)   LIST 

LIST  =     WEPE  NOT  hMJKJ^Y  TFWFPED.    IN  ORDER  TO' 
CALL  QCMDV(15,17) 
UFITE(6,805)   LIST 

LIST  =    'PRE\^NT  TTE  PROGRAM  FROM  CRASHING  VDU' 
CALL  QCMJw'(15,16) 
WRITE(6,805)    LIST 
LIST  =    'MUST  START  AGAIN.' 
CPU.  GOCV(15,15) 
WRITE(6,805)    LIST 

LIST  =     PFESS  ANY  KEY  TO  CONTINLE' 
CALL  QCMDW(15,13) 
WRITE(6,a05)   LIST 
CALL  QIM<EY(I3,I4) 
GOTO  20 
EM)IF 
C 

C  ALL  REQUIRTENTS  Hfl^G  ffiEN  SATISFIED.      LI^EAR  FEERESSION   IS 

C  PERFORMED  ON  TFE  TAFEL  FE6I0NS  AM)  TTE  RESULTS  ARE  (FJV5FVED 

C 

Dl  =    'Bc=    ' 
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YN«-e  =     V^DLTS' 

DO  35    I    =   1,01 

XDAT(I)    =  t«ATA(I,2) 
YDAT(I)    =  r«ATA(I,l) 
35  COVfTINUE 


X2  = 

0 

XI  = 

0 

XY  = 

0 

NI  = 

ARR(C2-C1)  +1 

YI  = 

0 

CZ   = 

J- 

-(4-k:i) 

C4  = 

J- 

-(4-K2) 

m  40  ] 

:  =  C4,C3 

XI 

= 

XI  ^-  r«ATA(I,2) 

X2 

= 

X2  +  (M)ATA(I,2)*t2) 

XY 

= 

XY  -•-  (MQATA(I,2)«MDATA(I,1)) 

YI 

= 

YI+  MDATA(I,1) 

40  CCNFIMJE 

^e  =  FUQAT(NI) 

Bl   =    ((XY»XI)-(YI«X2))/((XI«XI)-(^e«X2)) 

e2  =  ((XY«re)-(Yi*xi))/((X2»NS)-(xi»xi)) 

YI   =  M3ATA(J,1) 
XI   =    ((Y1-B1)/B2) 
X2  =  rf>ATA(C4,2) 
DX   =    (X2-Xl)/(50.) 
DO  45   I    =   1,50 

YDATl(I)    =    (B2<X2)+B1 

XDATl(I)    =  X2 

X2  =  X2-  DX 
45  CUvfTIMJE 

B2  =  -B2 
GZALL  CI^^/(D2,B2) 
X2  =  O 
XI    =  O 
XY  =  O 
NB  =  O 
YI    =  O 
A3  =  J+4-^Al 
A4  =  J+4+/C 
DO  50    I    =  A3,A4 

XI    =   XI    +  rf?.~TA(I,2) 

X2  =  X2  +    (M}ATA(I,2)*«2) 

XY  =   XY  +    (M)ATA(I,2)«M)ATA(I,1)) 

YI    =   YH-  M)ATA(I,1) 
50  CONTINUE 

NI   =  ABS(A3-A4)+1 

NS  =  n_QAT(NI) 

Bl   =    ((XY*XI)-(YI»X2))/((XI«XI)-(NB*X2)) 

B2  =    ((XY*NB)-(YI«XI))/((X2«NS)-(XI«XI)) 

YI   =  ro)ATA(J,l) 
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XI  =    ((Y1-B1)/B2) 
X2  =  M>ATA(A4,2) 
DX   =    (X2-Xl)/(50.) 
DO  55   I    =  51,100 

YDATl(I)    =    (B2«X1)^^1 

XDATld)    =  XI 

XI  =  X1+  DX 

55  CCNTINLE 

D3  =    '    Ba=    ' 

CflLL  CI1V^'(D4,B2) 

X2  =    (XDATl(50O)+XDATl(5Ol))/2 

TI71_E2=D1//D2//D3//D4 

^G  =   100 

CfiLL  GF¥«^l(XDAT,YDAT,01,Jl,^lfi^E,XDATl,YDATl,^e,Y^JPrE,TITl_E2 
+  ,A1,A2, CI, C2, TITLES) 

C 

C  THIS  PftFTT  PLLOWS  TVE  USER  TO  SELECT  EITTER  ANDD-ER  fNJDlC  CJR 

C  CATHDDIC  J^^FU.  REGION  OR  TO  REVIEW  Tl-E  GWFH  WI7HXJT  HftVTiVG 

C  TO  E6TER1I^e  BOTH  TVE  TfiFEL  REBIOhG  PGP\IN. 

C 

56  CPLL  QCLEftR(0,2) 
LIST  =    '     ' 

CPU.  QCMCV(20,16) 
li«ITE(6,805)    LIST 

LIST  =    •  IF  TVE  GRflPH  DID  MJT  SHOW  A  GOOD  CORREL/^TION' 
CALL  000^(20,15) 
WRITE(6,805)    LIST 

LIST  =      YOU  CAN  FEFEAT  TVE  CALOJLATIOre  FOR  TVE' 
CALL  aCMCV(20,14) 
WRITE(6,a05)   LIST 
LIST  =      BRATVCH   IN  QUESTION- 
CALL  QCM0V(20,13) 
WRITE(6,a05)   LIST 
LIST  =      1.   REPEAT  CALCULATIONS- 
CALL  QCM>/(20,11) 
WRITE(6,805)   LIST 
LIST  =    -2.   EXIT  THIS  PORTION- 
CALL  acrfDV(20,9) 
WRITE(6,805)   LISF 
LIST  =     ENTER  A  1   CF  A  2- 
CALl.  QCM^(20,7) 
WRITE(6,805)   LIST 
CALL  QCMDV(25,6) 
READ(6,aiO,EF^R=56)    M 
IF    (M    .ED.    1)    TVEN 

GOTO  30 
EM)IF 
IF    (M    .EQ.    2)    TVEN 

GOTO  20 

GOTO  56 
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EMJIF 
EM)IF 

C  TVE  CENTRAL  DIFFERMZE  I-E1H3)   IS   TO  BE  USED.    THE  CFTION  PROCEEDS 

C  IN  71-E  SPtE  MPlvrER  AS  OPTION  2 

IF  (M  .EQ.  3)  TVEN 
Al  =  O 
A2  =  O 
CI  =  O 
C2  =  O 
57  CALL  SLDFt(MDftTA,01,Jl,Nflft,Al,A2,Cl,C2,Sl,S2,S3,S4) 

IF    ((CI    .6E.   C2)    .OR.    (Al    -GE.   A2) )    Tl-EN 

C:ALL  QSM3DE(3) 

C:flLL  QCLEAR(0,4) 

LIST  =    '     ' 

C:ALL  QCMJ^(15,20) 

WRITE(6,805)   LIST 

LIST  =     YOU  DID  NOT  Pn  FTT  BOTH  AN  fNODIC 

CALL  (3CMDV(15,19) 

UFITE(6,805)    LIST 

LIST  =  AND  CA7HX)IC  BRANCH  OR  ELSE  ThE  POINTS- 
CALL  G3Dl>/(15,ia) 

WRITE(6,805)  LIST 

LIST  =  WERE  NOT  PROPERLY  TRPPF^D.    IN  ORDER  TD' 

CALL  GCJi:V(15,17) 

l^ITE(6,805)  LIST 

LIST  =     PREVENT  THE  PROGPyW  FTOI  O^ASHir^e  YOU' 

CALL  acro/(i5,i6) 

WRITE(6,a05)   LIST 

LIST  =     MUST  ST^TT  AGAIN.' 

CALL  (3CM3V(15,15) 

l*?ITE(6,805)   LIST 

LIST  =    'FfESS  ANY  KEY  TD  CONTINLE' 

CALL  QD1CW(15,13) 

li*?ITE(6,805)    LIST 

C:ALL  (3IM<EY(I3,I4) 

GOTO  20 
EM)IF 
Dl='Bc:=    • 
YNA^E  =    'VO.TS' 
DO  60   I   =   1,01 

XDAT(I)    =  M)ATA(I,2) 
YDAT(I)    =  rOATAdjl) 
60  (XNTIMJE 

X2  =  O 
XI  =  O 
XY  =  O 

NI   =  ABS(C2-C1)    +1 
YI   =  O 
C3  =  J-(2+Cl) 
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C4  =  J-(2-K:2) 
DO  70  I  =  C4,C3 

XI  =  XI  >  riDATA(I,2) 

X2  =  X2  +  (MDATA(I,2)««2) 

XY  =  XY  +  (t«ATA(I,2)«M3ATA(I,i)) 

YI  =  YI+  M)ATA(I,1) 
70        CONTINLE 

ise  =  FUDATCNI) 

Bl  =  ((XY«XI)-(YI»X2))/((XI«XI)-(^6*X2)) 

B2  =  ((XY«VG)-(YI«XI))/((X2<WB)-(XI«XI)) 

Yl  =  MDATA(J,l) 

XI  =  ((Y1-B1)/B2) 

X2  =  M3ATA(C4,2) 

DX  =  (X2-Xl)/(50.) 

DO  90  I  =  1,50 

YDATl(I)  =  (B2tX2)+Bl 

XDATl(I)  =  X2 

X2  =  X2-  DX 
80        CONTINUE 
B2  =  -82 
CALL  CaW(D2,B2) 
D3  ='  Ba=  ' 
X2  =  O 
XI  =  O 
XY  =  O 
NS  =  O 
YI  =  O 
A3  =  J+2+A1 
A4  =  J+2+A2 
DO  90  I  =  A3,A4 

XI  =  XI  +  MDATA(I,2) 

X2  =  X2  +  (M)ATA(I,2)*«2) 

XY  =  XY  +  (hDATA(I,2)«M)ATA(I,l)) 

YI  =  YI+  M)ATA(I,1) 
90        CDNTIMJE 

NX  =  ABS(A2-Al)-»-l 

NB  =  FLQAT(NI) 

Bl  =  ((XY»XI)-(YI«X2))/((XI»XI)-(NB*X2)) 

B2  =  ((XY«ve)-(YI«XI))/((X2«NB)-(XI»XI)) 

Yl  =  M>ATA(J,1) 

XI  =  ((Y1-B1)/B2) 

X2  =  rOATA(A4,2) 

DX  =  (X2-Xl)/(50. ) 

DO  lOO  I  =  51,100 

YDATl(I)  =  (B2»X1)+B1 

XDATl(I)  =  XI 

XI  =  X1+  DX 
100       CDNTINUE 

CALL  CarsW(D4,B2) 

T 1 7l£2=Dl/ /D2/ /D3/ /D4 

NB  =   100 
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CflLL  GRfiFVU(XDAT,YI>^T,01,Jl,JVWrt,XDATl,YDATl,re,YNPft,TITLE2 
+  ,A1,A2,C1,C2,TITLE3) 

58  CflLL  QC1J£AR(0,2) 

Lisrr  =  '    • 

CPLL  QCMQV(20,16) 

WRITE(6,805)    LIST 

LIST  =    '  IF  TVE  6F¥*=H  DID  NJT  SHDW  A  GDDD  CCFaFELATION' 

CPU.  aCMDv'(20,lb) 

WRITE(6,a05)    LIST 

LIST  =     YOU  CPN  I^PEAT   THE  CPLOJLATICWS  FOR  TTE' 

CPLL  QC^l]^/(20,14) 

WRI7E(6,805)   LIST 

LIST  =     BRfiNCM   IN  QLESTIGN' 

CPLL  aCMC]V(20,13) 

WRITE(6,805)   LIST 

LIST  =1.    FEPEAT  CPLOJLATIONS' 

CflLL  00*0^(20,11) 

l*?ITE(6,805)   LIST 

LIST  =    '2.    EXIT   THIS  PCFTIGN' 

CflLL  QD1DV(20,9) 

WRITE(6,805)    LIST 

LIST  =     ENTER  A  1   DR  A  2' 

CflLL  QCM>/(20,7) 

WRITE(6,805)    LIST 

CPLL  QCM>/(25,6) 

READ(6,810,ERR=5e)    M 

IF    (M    .EQ.    1)    TFEM 

GDTD  57 
EM)IF 
IF    (M   .ED.    2)    TVEN 

GOTO  20 
ELSE 

BDTD  58 
EJOIF 
EM)IF 

C  CPTION  4  RETLFre  THE  USER  TO  TVt  BEGI^^^I^G  CF  TVE  PRDGFVW 

IF    (M    -EQ.    4)    Tl-EN 

GOTO   10 
EM)  IF 

C  CPTICW   5   IS  TVE  O^ER_AY.    CN_Y  TFE  POTENT  I  ODYNPtllC  OJAES  CPN 

C  BE  QVRHED  GR  PLOTTED.    ThE  ORIGINflL  FILE   IS  STILL  AVAILABLE 

C  FOR  FIRTT-ER  «Y^JIRJLATIGN 

IF  (M  .EQ.  5)  THEN 
TrnjE3='  ' 
105       CALL  aa.EAR(0,5) 
LIST  = 
CALL  QCrtV(15,20) 
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URITE(6,B05)  LIST 

CPU-  QC11CW(25,12) 

LIST  =  ENTER  Tl-E  FILE  Nf^TE  TD  BE  U^D' 

WRITE(6,a05)  LIST 

CALL  QCMjM(20,11) 

LIST  =     A  BLftK  TD  EXIT,    CR   "DIR"   FDR  A  DIRECTORY' 

l4RITE(6,a05)   LIST 

CALL  QCMDV(36,10) 

READ(6,80O)   NP^E2 

CALL  aCMDV(40,9) 

IF    (NPTEZ    -EQ.    •  '  )    T^€N 

QQ7D  21 
END  IF 

IF    (NPrE2    .ED.      DIR)    TVEN 
CALL  aCL£AR(0,5) 
CALL  aCMDV(l,23) 

PAUSE     ENTER   "DIR   *.    /W"   OR  A  BL/¥^  TO  CONTIME' 
CALL  aLECRL(23,6, 1,1,79,0,5) 
GOTO   105 
EM)  IF 
DO  750   I   =  1,01 

XDAT(I)    =  «)ATA(I,2) 
YDAT(I)    =  M]ATA(I,1) 
750  ODNTIMJE 

CALL  DATAIN(CDATAl,02,NfifC2,IER) 
IF    (lER    -ED.    1)    TI-EN 

GOTO   105 
EJVOIF 

CALL  DATAM(CDATAl,MDATAl,a2,J2) 
CALL  a-En<(n>ATAl,a2,J2,Sl,S2,S3,S4) 
DO  760   1=1,02 

XDATl(I)    =  rt)ATAl(I,2) 
YDATl(I)   =  M)ATA1(I,1) 
760  CDNTIMJE 

TITLE2  ='(X€R_AY      //NfiTE//*     ' //Nflrt2 
YNATE  =     VOLTS' 
NAr€2=NAft 
NPtE  =    '     ' 

CALL  GRAPHl(XDAT,YDAT,01,Jl,Nflrt,XDATl,YDATl,a2,YNArE,TITL£2, 
+  A1,A2,C1,C2,TITLE3) 

NPliE  =  Nfift2 
CALL  aCLEAR(0,2) 
GOTO  21 
ENDIF 

IF    (M   .EQ.    6)    Tl-EN 

GOTO  900 
ELSE 

GOTO  20 
ENDIF 
80O  FORMAT(Ae) 

805  FCFJ1AT(A50,\) 
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810  FDRMATdl) 

900  Cftl.  GEI1DDE(3) 

EM) 
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»        THIS  SLjBFOJTI^E  fEPDS  TVE  DATA  FILE  6BVERTATED  BY   Tl-E 
«        MODEL   351  UHICH  HAS  BEEW  TFiPiseFEFiRED  OnO  A  DISK 

C  CDATA  =  TVC  AFRAY  WHICH  CC^^■AINI^E  PG^^E^fTIAL  (=^0 

C  OJRFBNfT  DE^6ITY 

C  NPIt  =  TVE  rv^Vt  GF  TVe  DATA  FILE 

C  CI  =  CHARACTER  STRI^e  RESENT    IN  Tl-E  DATA  FILE 

C  C2  =  CHAFWCTER  STRI^e  KtaJVfT    IN  TFE  DATA  FILE 

C  N  =  NUMBER  CF  DAFA  POINTS   IN  Tl-E  FILE, 

C  DETER11^ED  DtJRIhG  T>E  PRDGRf^ 

SUBROUTINE  DATAIN   ( CDATA, N, NATE,  lER) 

HEAL«8  CDATA(2000,2) 

INrEGER»2,    N,I 

OHP^ACTTER  C1»17,C2<6,NP^E«8,LIST*50,C3«50 

OPEN   (UNIT  =   lO,   ERR  =  200, FILE  =  NATE,   STATUS  =     OLD'  ) 

DO  50  N=l,2000 

READ  (10,900,EFR=60)  C1,CDATA(N,1)  ,C2,CDATA(N,2) 
50     CDNTINLE 
60     FEWIM)(10) 

DO  70  I  =  1,2000 

READ(10,«,ERR  =  100) 
70     CONTIMJE 
100    IF  (I  .NE.  N)  TFEN 

CALL  QSMDDE(3) 

CALL  QCLEAR(0,4) 

LIST  =  '  ' 

CALL  QCrf^/(15,20) 

WRITE(6,905)  LIST 

LIST  =  AN  ERROR  HAS  BEEN  EM3XWTERED  IN* 

CALL  QCrfDV(15,19) 

WRITE(6,905)  LIST 

LIST  =     THE  FCFWAT  OF  TFE  DATA  FILE.' 

CALL  QCrDV(15,18) 

WRITE(6,905)    LIST 

LIST  =  THIS  LEUPLLY  ir**a_V^S  A  SWITCH  BETVEEN' 

CALL  QaiOV(15,17) 

WRITE(6,905)  LIST 

LIST  =     FLOATING  HDINT  AND  EXFCTENTIAL  MJTATION' 

CALL  QCMCW(15,16) 

WRITE(6,*K)5)    LIST 

LIST  =  AROUM)  ECOOR. ' 

CALL  001*0/(15,15) 

U«ITE(6,905)  LIST 

LIST  =  TO  CORFECT  THIS  FRDBLEM  EXIT  TFE  PROGFV^I' 

CALL  aCM0V(15,13) 

WRITE(6,905)  LIST 

LIST  =  AM)  USE  TTE  DOS  "EDLIN"  FEATLFE.  TVE' 
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Q^LL  aCMDW(15,12) 

WRITE(6,905)  LIST 

LIST  =  PRDELEM  aEBi^G  AT  LINE   #: ' 

CALL  QDiL}V(15,li) 

URIfE(6,905)  LIST 

CPLL  QCMJ/(38,11) 

U«I7E(A,915)  N 

LIST  =     PFESS  PNY  KEY  TD  CCI4TINLJE' 

CPLL  GD-C]V(15,8) 

WRITE(6,905)  LIST 

CPLL  QIM<EY(I3,I4) 

GOTO  990 
EM)IF 
GOTO  990 
200  CflLl_  aSMXE(3) 

CflLL  Qa-£PR(0,4) 
LIST  =  •  ' 
a:^LL  aCMDk"l'-,20) 
WRI7E(6,905)  LIST 

LIST  =    'TVE  FILE  NfiTE  ENTEFED  DOES  NOT  EXIST.' 
CPLL  aCM3V(15,19) 
WRI7E(6,905)   LIST 

LIST  =     TO  EXPmi^  TT-E  DIFEETCRY  OF  POSSIBLE  FILES' 
CflLL  QC11DV(15,18) 
WRITE(6,905)    LIST 

LIST  =    'ETvrrER   "DIR"   W-EN  ASKED  TO  ENTER  ThE  FILE' 
CALL  QCMDV(15,17) 
WRI7E(6,905)   LIST 

LIST  =    'IS¥»t  FOLLOWED  DY  TFE  DOS  CCftlAM). ' 
CALL  QCMJV(6,905) 
U«ITE(6,905)    LIST 
LIST  =    '"DIR  «.    /W" 
CALi.  QCM]V(35,16) 
WRnE(6,905)   LIST 

LIST  =      PRESS  ANY  KEY  TO  COrTTIMJE' 
CALL  QCM)V(15,14) 
WRITE    (6,905)   LIST 
lER  =  1 

CALL  QIM<EY(I3,I4) 
GOTO  990 
900  FOFWAT    (A17,GB.4,A7,610.4) 

905  FDFMAT(A50,\) 

915  FCR1AT(I4,\) 

990  RERJRM 

END 
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«      THIS  SUBFaJTI^E  MfiNIRJLATES  Tl-E  DATA   IN    TVE  PffJAY  CDATA 
«      T»-E  PRCFER  FORI  AM)   IDTENIFIES  EOTJR 

C  CDATA  =   THE  f^R^Y  CC^^■AINI^G  POTHNfTAIL  ATO  CURRENT 

C  DBNBITY 

C  rtJATA  =  CDATA  CDT^/ERTED  TO  POTENTIAL  fM}   LOG 

C  CURRENT  DENSITY 

C  N  =  NLTBER  OF  DATA  POINTS  IN  TTE  DATA  FILE 

C  J  =  LOCATION  WITHIN  r«)ATA  WHERE  ECORR  IS  LOCATED 

C  I  =  DO  LOOP  COLNTER 

C  XI  =  ABSOLUTE  \/AUJE  OF  CURRENT  DENSITY  AT  I 

C  USED  FOR  DETEFMINING  EC^R 

C  X2  =  ABSOLUTE  WLUE  OF  CURRENT  DENSITY  AT  I-i 

C  USED  FOR  DETERMINING  ECORR 

C  X3  =  LOGiO(Xi) 

^UEFOJTI^E  DATAM    (CDATA,rCATA,N,J) 
REAL*8  MDATA(2000,2),CE)ATA(2000,2) 
INTEGER«2  N,J 
X2  =   100000. 

C  OJhMERT  CLF«ENT  DENSITY    INTO  LOB  CURRENT   DENSITY 

DO  100   I=l,N-2 

XI  =  ^«S(CDATA(I,2)) 
X3  =  LOBIO(XI) 

C      IF/TVEN  TEST  UBED  IN  DETERMINING  LOCATICN  OF  ECOOR 

IF  (  XI  .LT.  X2)  TVEN 
J  =  I 
X2  =  XI 
ELSE 

X2  =  X2 
ENDIF 

M)ATA(I,1)    =  CDATA(I,1) 
M)ATA(I,2)    =   X3 
100  CONTINUE 

RETfLHM 
EM) 
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THIS  SLJBHaJ^I^e  FtFFOFWS  7VO  FUNCTIOve.    ITS  FIRST  FUNCTIGN 
IS  TO  OECK  TO  EN3UPE  1V¥^J  Tl-E  POTENTIPL  VALLJES  RECDRDED 
ACOJRftTELY  REFLECT  TKE  rtmjflL  POLPRIZATIGN.    TVE  SECOND   IS 
TO  r}ETE^WI^E  TVE  MJTBER  of  POI^fTS   IN  THE  PNUDIC  AM)  CATHDDIC 
BRflM>ES  WHICH  MftY  HE  USED   IN  TVE  CENTFifiL  DIFFEFNCE  fETHDD. 

rCATA    :    THAT  AFRAY  COVfTAINir^B  POrTE^^"IAL  AW)  LOG  CURRENT 

DE>eiTY 
N3  :    TVE  NUTBER  GF  POINTS   IN  TVE  CATVDDIC  BRPlvCH  WHICH 

CAN  BE  USED   IN  TVE  CE^fTRflL  DIFFERENCE  rEHVCD 
N4  :    TVE  NLTBER  OF  POINTS   IN  TVE  PNODIC  BRP^CH  WHICH 

CAN  BE  USED   IN  TVE  CENTTV»_  DIFFERCE  rtTVDD 

OTHER  V/PRIBLES  ARE  USED  rERELY   IN  CVECKING  REQUIFtD  CDNDITICIVB 

SUBRDUTI^E  CVEa<(M)ATA,N,J,Nl,N2,N3,N4) 
REAL«8  M)ATA(2000,2) 
INTEGER«2   I  ,J,N,N1,N2,N3,N4 
CHfiRACTER«50  LIST 

C  TVE  FIRST    INCRErevfT   IN  POTENTIAL  FOR  TVE  CATHDDIC  BRAISCH      IS 

C  DETERMIhED 

lO  Y  =  MDATACJ,!) 

Yl   =  rC)ATA(J-l,l) 

Zl  =  ABS(ABS(Y1)-ABS(Y)) 

Zl  =  ANINT(Z 1/0. 0001) 

C  SUBEEGUENT  POIENTIAL   INCREI'ENTS,    WHICH  MJST  BE  EOmL  TO  TIE 

C  FIRST  ARE  CVECKED.      LOG  dJVJREMT  CENBITY  W^ULES  AFE  PLSD  CVECKED 

C  TO  ENSURE  THAT  TVEY  WILL  NOT  RESULT   IN  TVE  CALCULATION  OF  A 

C  SLGFE  OF    I^FINITE  \/PLUE. 

DO  20   I    =  J, 2,-1 
Y  =  MDATA(I,1) 
Yl   =  MDATA( 1-1,1) 

X  =  M)ATA(I,2) 

XI  =  MDATA( 1-1,2) 

Z  =  ABS(ABS(Y1)-ABS(Y)) 
Z  =  ANINT(Z/0.0001) 

C  OTE  CF  TVE  TVD  FJEQUIRED  COMDITIOTse  VWS  FAILED. 

IF    ((Z    .^E.    Zl)    .OR.    (    XI    .LE.    X    ))   GOTO  25 

20  CONTIMJE 

25  Nl   =    1+2 
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C  TJ-E  FIRST    INCREJ-EMT    IN  POTtNTI<=4-  FOR  TKC  fitODlC  BRfiNDH    IS 

C  iDETER1I^eD 

X   =  f«)ATA(J+l,l) 

Y  =  r"«)ATA(J+2,l) 
Zl   =  ABS(X-Y) 

zi  =  <^ir^n"(zi/o.oooi) 

C  SLBSEOJENT  POTENTIAL   INCFOtNrS,    UHICH  MUBT  BE  EQUAL  TO  THE 

C  FIRST  Pf^  a-EC3<ED.      LCG  OJFRENT  D6^BITY  VALUES  ARE  <3LS0  OECKED 

C  TO  ENSURE  T>*^T  T>EY  WILL  NOT  RESULT    IN  THE  CALCULATION  CT  A 

C  SLOPE  OF    I^FINITE  VfiLUE- 

DO  30   I   =  J+2,N-1 

Y  =  |-t)ATA(I,l) 

Yl   =  rfDATAC  1+1,1) 

X  =  rCATA(I,2) 

XI  =  M)ATA(  1+1,2) 
Z  =  ABS(Yl-Y) 

Z  =  ANINT(Z/0.0001) 

C  OVE  OF  TVE  TWO  REQUIRED  COM)ITICfve  HAS  FAILED. 

IF    ((Z    .^E.    Zl)    .OR.    (    XI    .LE.    X    ))    BOTO  35 

30  CONTINLE 

35  N2  =    I    -  2 

N3  =  J   -  Nl 

N4  =  N2  -  J 

C  TFE  NUnSER  OF  POINTS  Pt^lLfOLE  FCF  EACH  BRflNCH   IS  COTPPRED 

C  WITH  A  MIMJIIUM  NUTEER  OF   15.    IF  ^EITVER  BRflM»  HfiG  ENOUGH 

C  AVAILABLE  POINTS  VARIOUS  rESSAGES  ARE  DISPLAYED  LISTING 

C  POSSIBLE  CALJBES  PIM)  CFTIONS. 

IF    ((N3    -LE.    15)    .OR.    (N4    .LE.    15))    TVEN 
CALL  QCL£PR(0,4) 
LIST  =    •     ' 
CALL  QCMCV(20,20) 
WRITE(6,805)    LIST 

LIST  =      TVE  CENTFWL  DIFFERVCE  rEThtD  REQUIRES  THAT' 
CALL  QaiDV(20,19) 
URITE(6,805)   LIST 

LIST  =     TVE  CURRENT  BE  rEABJFED  AT  EQUf^L  DIFFERENTIALS' 
CALL  QD1CV(20,18) 
WRITE(6,a05)   LIST 

LIST  =     OF  TVE  AFT^IED  VOLTAGE.    TVE  RfiR  MDEL  351' 
CflLL  QCrt:V(20,17) 
WRITE(6,805)    LIST 

LIST  =     RECOFBS  TVE  CURRENT  AT   IfCREJIENTS  OF  2.0  nA/' 
CALL  aCM>/(20,16) 
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URITE(6,B05)    LIST 

LIST  =     UNTIL  T>E  CHftM3E   IN  OFRENT  EXCEEDS  A  PRESCRIBED' 
CnLL  QCM>^(20,15) 
WRITE(6,a05)    LIST 

LIST  =     LEVEL.    AT  TT^T  POINT  Tl-E  351   BEGI^e  RECQRDI^G• 
CALL  QCI-CV(20,14) 
l«RITE(6,a05)    LIST 

LISrr  =     TVE  CURRENT  AT   INCFQ-evfTS  OF  0.5  oA/.   AT  AN' 
CPUL  QCM3V(20,13) 
WRITE(6,805)    LIST 

LIST  =    'HIQ-ER  APPLIED  Va_TAGE  7TE  CHPISCE   IN  CURRENT' 
CALL  QaiGV(20,i2) 
WRITE(6,805)    LIST 

LIST  =    'DECREASES  TO  A  POINT  WHERE  TTE  CURRENT    IS  PCAIN' 
CALL  QCM7v'(20,ll) 
WRnE(6,a05)    LIST 

LIST  =    '^E/::«JRED  AT   INCRETtNIS  OF  71-E  PFR_IED  VOLTAGE' 
CALL  QCM3V(20,10) 
WRITE(6,805)    LIST 
LIST  =     EGUfiL  TO  2.0  ntV' 
CALL  aCM3/(20,9) 
WRITE(6,805)    LIST 

LIST  =     PRESS  ANY  KEY  TO  CONFINUE' 
CALL  aa«V(20,7) 
WRITE(6,805)    LIST 
CALL  QIIM<EY(I,r«) 
CALL  QCLEAR(0,4) 
400  LIST  =    '     ' 

CALL  aCMDV(20,20) 

V«nE(6,805)    LIST 

LIST  =  Tl-E  DATA  RECORDED  FOR  THIS  EXPERirtNT  EITVER: ' 

CALL  QCMDV(20,19) 

WRITE(6,a05)  LIST 

LIST  =  '  ' 

CALL  aCM3-'(6,18) 

WRITE(6,805)  LIST 

LIST  =  '1.  FCFLECTS  THAT  Tl-E  CHANGES  IN  CLFRENT' 

CALL  QCMDV(20,17) 

URITE(6,805)  LIST 

LIST  =  'MEAeUFE^ENT  DESCRIBED  EARLIER  DID  NOT  OCCUR' 

CALL  acrOV(20,16) 

URITE(6,B05)  LIST 

LIST  =  AT  A  POINT(S),  F^LATIV^  TO  ECODR,  THAT  WOULD' 

CALL  QaiDV(20,15) 

WRITE(6,a05)  LIST 

LIST  =  'GE^eRATE  ENOUGH  DATA  POINTS  FOR  TFE  USE  OF  TVE' 

CALL  QCM^(20,14) 

l««RITE(6,805)  LIST 

LIST  =     CENTTVl.  DIl-hhhtNCE  l*ET>CD. ' 

CALL  QCrtV(20,13) 

URITE(6,S05)   LIST 
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LIST  =  OR- 

CPLL  QCrOV(39,ll) 

WRH1E(6,805)  LIST 

LIST  =    '2-    TVE  MTOEL  351   DID  NOT  FtCCF?D  TT-E  APPLIED' 

CPLL  QD"tJw'(20,9) 

WRITE(6,805)    LIST 

LIST  =  POTENTIAL  PRDPERLY.  THIS  CP>N   BE  CORFECTED' 

CALL  00*0^(20,8) 

WRnE(6,805)  LIST 

LIST  =    '  IF  TTC  PROBLEM   IS  NOT  COFFECTED  TTEN  TVE' 

CflLL  QCM]V(20,7) 

U«ITE(6,805)   LIST 

LIST  =     CENTRA  DIFFERNCE  fETHDD  WILL  NOT  UCR<,    AM)  THE' 

CALL  aCMDV(20,6) 

WRnE(6,805)   LIST 

LIST  =     CUBIC  SR_I^E  fETHOD  MAY  NOT  UCR<' 

CALL  QCMDV(20,5) 

wRrrE(6,ao5)  list 

LIST  =  press  any  key  TO  CCNTINUE' 

CALL  00*0^(20,3) 

URITE(6,S05)  LIST 

CALL  OIM<EY(I,MB) 

CALL  QCLEAR(0,4) 

LIST  =    '  IF  AN  ATra-PT   IS  MADE  TO  CORFECT  TKE' 

CALL  QOC1V(20,17) 

V«FITE(6,a05)    LIST 

LIST  =     PRCBLEM  JhE.  CORRECT  fETHD   IS  DETER1INED' 

CALL  ac^c^'(20,l6) 

WRnE(6,a05)    LIST 

LIST  =    '  IF  THE  ATTEfFT    IS  UTBLCCESSFUL  TFEN  Tl-E' 

CALL  QCMDV(20,15) 

WRITE(6,805)    LIST 

LIST  =     CUBIC  Sn_IhE  rtTHDD  SHOULD  BE  USED.' 

CPLL  QCMDV(20,14) 

WRITE(6,805)    LIST 

LIST  =    '  IF  THE  ATTE^PT   IS  SUCCESSFUL  FITTER' 

CALL  QCrCV(20,13) 

l»RITE(6,B05)    LIST 

LIST  =     rETHOD  MAY  BE  USED.' 

C:ALL  GCMDV(20,L2) 

l««ITE(6,EX)5)    LIST 

list  =  '  if  this  set  of  ressages  appears  ag^in' 
c:all  qcmdv(20,io) 
write(6,805)  list 
list  =    after  ariem^ing  to  correct  tl-e  problem' 

CALL  aCM0V(20,S') 
WRITE(6,805)  LIST 
LIST  =  TFEN  TFE  CLFRENT  DENSITY  V«_UES  ff^   SUCH' 

CALL  QCM0V(20,8) 
li«ITE(6,a05)  LIST 
LIST  =  THAT  ^E1T^ER  rETHDD  CfW  BE  USED' 
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CPU-  aCMDV(20,7) 
WRITE(6,805)    LIST 

LIST  =      PRESS  PNY  KEY  TO  CGNTIMJE' 
CALL  a:MDV(20,5) 
WRITE(6,a05)    LIST 
CflLL  QIM<EY(I,MB) 
420  CflLL  QCLJD=«(0,2) 

LIST  =      1.    DO  MDT  CORRECT  Tl-E  PRDBLEM. ' 
CPLL  aCMDV(25,13) 
WRITE(6,B05)    LIST 

LIST  =    '2.    "mv  PMl>  correct  TVE  PRDBLEM,' 
CPLL  aCM3V(25,ll) 
WRITE(6,a05)   LIST 
LIST  =     ENTER  A   i  CF?  A  2' 
CALL  QCMDV(32,9) 
WRITE(6,805)   LIST 
CALL  aCM0V(39,8) 

REA0(6,^'.  ,ERR  =  450,    lOSTAT  =  J4)    M 
450  IF    ((M    .GT.    2)    -OR.    (J4    .^e.    O))    Tl-EN 

CALL  aCL£AR(0,4) 

LIST  =    '     ' 

Cfll_  acrtDV(21,16) 

WRITE(6,805)    LIST 

LIST  =     YOU  DID  ISDT  ENTER  AN   INTEEER  OF' 

CALL  QCMDV(20,15) 

WRrTE(6,805)    LIST 

LIST  =     V/PLLE   i  OR  2' 

CALL  aCI1DV(37,14) 

WRITE(6,905)   LIST 

LIST  =  RRESS  ANY  KEY  TO  CONTINUE' 

CAUL  aCMDV(25,12) 

URITE(6,805)  LIST 

CALL  aCMDV(39,ll) 

CALL  QIM<EY(I,rB) 

GOTO  420 
EM)IF 
IF    (M   .ED.    I)    TVGM 

N9  =   1 

GOTO   lOOO 
ETOIF 

C  AN  ATTEr-PT    IS  BEING  MADE  TO  CORRECT  TVE  POTEf^AIL  VALUES  SUCH 

C  THAT  TT-EY  ACCURATELY  REFLECT   ThE  EXPERIMENTAL  POTENTIALS. 

C  ELOCKS   500  AND  600  CORFECT   TTE  PROBLEM,    BUT  Gr4_Y  QTE  OF   TFE 

C  TiOD  BLOCKS  ARE  USED  FOR  ANY  Gl^JEN  FILE. 

IF    (M    .EQ.    2)    Tl-EN 
DO  500    I    =   1,N,2 
Yl   =  MDATA(I,1) 
Y2  =  M)ATA(  1+1,1) 
IF    (Yl    .ED.    Y2)    TT-EN 
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MDATA(I,1)    =  MDATAd,!)    -  O.00O5DO 
N9  =    L 
QViDIF 
50O  CONTIMJE 

DO  600   I   =  2,N,2 
Yl  =  M)ATA(I,1) 
Y2  =  M)ATA(  1+1,1) 
IF    (Yl    .ED.    Y2)    TVEN 

M3ATA(I,1)    =  M>ATA(I,1)    -  0.0005D0 
|Vf?  =   1 
EIOIF 
600  CDNTIMJE 

IF    (N9    .ED.    1)    TVEN 
CALL  aCLEflR(0,2) 
LIST  =    •     • 
CPLL  00-0/(30,14) 
WRITE(6,805)    LIST 
LIST  =     CLFWdLIIOvl  SUCCESSFUL' 
CALL  acrfDV(30,13) 
WRITE(6,805)    LIST 
LIST  =     USE  EITHER  rtTVCD' 
CALL  aCMDV(30,ll) 
WRITE(6,805)    LIST 

LIST  =     PRESS  ANY  KEY  TO  CONfTINUE' 
CflLL  CCM]V(30,9) 
WRITE(6,805)    LIST 
CALL  QIM<EY(I,I-B) 
GOTO   10 
ELSE 

CALL  GCLEAR(0,4) 

LIST  =     CORRECTICW  UNSUCCESSFUL' 
CALL  QCriDV(30,13) 
WRITE(6,805)   LIST 

LIST  =     USE  EITVER  CUBIC  g^irC  METHX)' 
CPLL  QCMaV/(30,li) 
URITE(6,805)   LIST 

LIST  =     PRESS  ANY  KEY  TD  (XNTIMJE' 
CALL  (3Crt]V(30,9) 
WRITE(6,a05)    LIST 
CflLL  QIM<EY(I,MB) 
GOTO   lOOO 
ETvOIF 
EM)IF 
END  IF 
905        FDR1AT(A50,\) 
810        FCRIAT    (II) 
lOOO     END 
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XDAT 

YDAT 

XDATl 

YDATl 

TX 


TY 

TXl 

TYl 


7>t   X   V«_LJES  OF  1>E  FIRST  3Rff^ 

Tl-E  Y  V<^UE5  GF   1>E  FIRST  GJV5PH 

TFE   X  W^UEB  CF  TFE  3EDCND  GRfiFH 

"n-E  Y  WLUES  OF  TFE  SECGND  GRfiFH 

XDAT  AS  CF?IGINAL  PASSED  PM)  SCPLED.    UEED  TO  PLUDW 

ITE  USER  TO  RESCPLE  THE  AXIS  ftC  MAINTAIN  ALL  TTE 

Vf«_LJE5  CRIGINflLLY  PASSED 

USED  AS  TX  BUT  FOR  YDAT 

UEED  AS  TX   HJT  FCR   XDATl 

UEED  AS  TX   aJT  FCR  YDATl 


C  IN  THIS  SLJBROLJTI^E  PLL  FRDT-RTS  PfE  REPLIED  TO  BY  TTE  UEE  OF 

C  StFT  KEYS.    THIS  MEP#v6  THAT  A  CPR^IAGE  RETURVi    IS  NOT  USED. 

SUEROUTirC  SVFHK  XDAT,  YDAT,  N3,  J,  NATt,  XDAT'  ,^'OATl  ,N4,YNArt  ,T1TLE2, 
+A1 ,A2,C1 ,C2,TITLE3) 
REPL»4   XDAT(200O),YDAT(200O),XDATl(200O),YDATl(2O0O),F,Fl,F2 
REAL«4   TX(2000),TY(2000),TX1(2000),TY1(2000),X2(1),Y2(1) 
INTEGER*2  K2,K3,K5,K6,L,XL,XR,YL,YR,rB,Al  ,A2,C1,C2,M,  IA(80) 
INTEGER«2    I  ,J,N,K7,K8,K9,    JCCL(IO)  ,JRDW(10)  ,Nl,he,N3,N4,N5,N6,N7 
IN1EGER«2  L1,L2,L3,L4 

CHARACTER  NP^e«8,YNfi^E«5,FACT*6,XNPIt*12,TITLE»15,TITLE2«32 
CHPfifi[JTEFW28,    OJES,    TITLE3«50 
N7  =   1 
M  =  O 

N  =  MAX(N3,N4) 
N3  =  N3  -  3 
N4  =  N4  -  3 
IF    (N4    .LE.    O)    N4  =   1 

C  Tl-E  ARF^YS  ARE  SCALED  SUCH  THAT  Tl-E  Y  AXIS  V;«_IJES  ARE  DISCERNIBLE 

C  FR3-I  EACH  OD-ER 

F  =  MIN(ABS(YDAT(1)),ABS(YDAT(N3))) 
Fl   =  mN(^«S(YDATl(l)),ABS(YDATl(N4))) 
IF    (N4    .GT.    2)    TTEN 

F2  =  MIN(F1,F) 
ELSE 

F2  =  F 
EJSOIF 
IF    (F2   .6E.    1.)    TVEN 

F2  =   1.0 

FACT  =      x   1' 
EM)IF 
IF    ((F2    .BE.    .1)    .PtM).    (F2    -LT.    1.))    TVEN 

F2  =   lO 

FACT  =      x    10' 
E^DIF 
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IF    ((F2    .GE.    .01)    .PND.    (F2    -LF.    .1))    TVEN 
F2  =   100. 
FACT  =      x    100' 
ENDIF 

IF    ((F2    .GE.    .001)    .PM).    (F2    .LT.    .01))    TI-EN 
F2  =  1000 
FACT  =      x   1000' 
ELSE 

IF    (F2    .LE.    0.0010)    Tl-EN 
F2  =   1000 
FPCT  =      x    1000' 
EM)IF 
ENDIF 
DO  6   I   =  1,2000 

YDAT(I)    =  F2«YDAT(I) 
YDATl(I)    =  F2*YDAT1(I) 
6  CXivfTIhLE 

5  N5  =  N3 

N6  =  N4 

C  71-E  PERMANENT  flFRAYS  PfE  FILLED 

DO  10  I   =  1,2000 
TX(I)    =  XDAT(I) 
TY(I)   =  YDAT(I) 
TXl(I)    =  XDATl(I) 
TYl(I)    =  YDATl(I) 

10       ajvrriNUE 

C  Tl-E  a^flPH  AXIS  AND  Mlt^XM  ARE  DETEFWI^ED 

XMAX  =  -1000 

XMIN  =   1000 

YMAX  =   XMAX 

YMIN  =   XMIN 

XMAXl   =  XMAX 

XMINl   =  XMIN 

YMAXl   =   XMAX 

YMINl   =   XMIN 

JCQLl   =  90 

JCI1J2  =  610 

JRCWl   =  35 

JRDW2  =  300 

DO  200   I   =   1,N3 

YMAX  =  MAX(YDAT(I),YMAX) 

XMAX   =  MAX(XDAT(I),XMAX) 

XMIN  =  MIN(XDAT(I),XMIN) 

YMIN  =  MIN(YDAT(I),YMIN) 
200        CXNTINUE 
205        DO  210   I   =   1,N4 

XMAXl  =  MAX(XDAT1(I),XMAX1) 
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YMAXl  = 

=  MAX(YDAT1(I),YMAX1) 

XMINl  -■ 

=  MIN(XEAT1(1),XMIN1) 

vniNi  •- 

=  MIN(YDAT1(I),YMIN1) 

210 

CONTIMJE 

215 

YFIN  = 

(MAX(YMAX,YMAX1))  + 

.1 

YSI   = 

(MIN(YMIN,YMirvU))  - 

.1 

XST   = 

(MIN(XMIN,XMI^a))  -. 

1 

XFIN  = 

(MftX(XMftX,XMPO(l))  + 

.1 

225 

YMIN  = 

YST 

YMAX  = 

YFIN 

XMAX  = 

XFIN 

XMIN  = 

XST 

YUMd  = 

YMIN 

XLKb  = 

XMIN 

XNOrt  = 

=  LUb  I  A/cm2' 

TIT!  F  = 

=  'PDTE^fTIaDYNftMIC• 

ICPT  = 

0 

YO^tHX 

=   0.3 

ASJ-tCT 

=  1.5 

TVE  GRPPH  AXIS  ARE  DRPWM  ftC  LABELED.    TITLES  PFfE  PLACED 

CALL  QSMjDE(16) 

CALL  QFLDT(JCa_l,JCaj2,JRDWl,JRDW2, XMIN, XMAX, YMIN, YMAX,XCR3,YDRB, 
+ICFT  ,YrPJERX ,  ASPECT ) 

XrV\IOR  =    (XF1N-XST)/10. 

XMIfvCR  =  O 

LABEL  =  -1 

M}EC     =  2 

YMPJGR  =(YFIN-YST)/10. 

YMIMJ?  =  O 

M<NT  =  5 

MXJT5  =  O 

IKOLDR  =  1 

ISY^BDL=-2 

KLJ^SYM=1 

CALL  (XIl_0R(O,O) 

CALL  QSETUP(^IXrrS,IKDLDR,ISY^BDL,Kl-RSYM) 

CPt-L  QXAXIS(XST,XFIN,XMfVJCR,XMINLJR,LABEL,NDEC) 

CPLL  OYAXISCYST, YFIN, YMPJCF,YMINCR, LABEL, NDEC) 

CALL  QXTICS(M<Nr,JCa_,vlROW) 

CAL:_  GGTXT(12,XNPrE, 2,340,1,0) 

CALL  aYTICS(|NMSfT,JCOL,JRGW) 

CALL  aGTXT(5,YNfrE, 2, 1,24^,-1) 

CALL  GGTXT(6, FACT, 2, 1,162,-1) 

CALL  CETXT(  15, TITLE, 2, 254, 335,0) 

CALL  GETTXTCS, NATE, 2, 382, 335,0) 

CALL  QGTXT(32,TITLE2,2,254,321,0) 

CALL  QGTXT(50,TITLE3, 2, 182,307,0) 

IKOLOR  =  2 

CALL  QSETUP(MXTT5,  IKOLOR,  I5YMBGL,KLRSYM) 
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CflLL  GRTDI(XST,YST,Li,L2) 
CPt±.  lF7DI(XFIN,YFIN,L3,L4) 
CjPLL.  C1_IM:(L1,L2,L1,L4,2) 
CPLL  GLIIVE(L1,L4,L3,L4,2) 
CALL  GLirvE(L3,L4,L3,L2,2) 
CfiLL  GI_I^E:(L3,L2,Li,L2,2) 

C  TVe  FIRST  CURVE    IS  DRfiUN.    Tl-E  LlhE  ai_GR   IS  CHANGED  PND  TVE 

C  SECCM3  CUF^E   IS  DR«WN 

CPLL  CTTAa_(l,N5,XDAT,YDAT) 

IKOLCF  =   1 

CflLL  QSETlJP(^OJTS,IKCL0R,ISY^1BDL,i<LJSYr1) 

CPU.  QTAa_(l,N£i,XDATl,yDATl) 

CPLL  rjBEEP 

M  =  M  +   1 

C  IF  TVE  CURVES  ARE  OF  DERIVATIVES  Pl^  A  LI^EAR  REGION  HAS  <=1J^EADY 

C  BEEN  SELECTED,    TVE  USER  MAY  DECIDE  TO  GE^ERATE  TVE  TAFEL  CONSTANTS 

C  BASED  ON  THIS  FOTTION  OF  TVE  OJRV/E. 

380  IF    ((    M    .ST.    1)    .AM).    (N4    .GE.O) )    TVEN 

QLES  =     GEIVERATE  TAFEL  CONSTANT' 

CALL  QGTXT( 28, QUES, 3,1,335,0) 

GUES  =  '      Y  OR  N  ?' 

CALL  QBTXT (28, QUES, 3, 1,320,0) 

CALL  QBEEP 

CALL  QIM<EY(I,M8) 

IF  ((l-B  .EQ.  89)  .OR.  (MB  .EQ.  121))  TVEN 

QUES  =  '  ' 

CALL  aGTXT(28, QUES, 3, 1,335,0) 

CALL  QGTXT(28, QUES, 3, 1,320,0) 

381  QUES  =  '  IS  THIS  TVE  ANODIC  OR' 
CAUL  QGTXT(28, QUES, 3, 1,335,0) 
CUES  =  CATVCDIC  BRANCH' 

CALL  aGTXT(28,aLES, 3, 1,320,0) 

OJES  =  •      A  OR  C  ?' 

C:ALL  QGTXT(28, QUES, 3, 1,305,0) 

CALL  QBEEP 

CALL  QIM<EY(I,MB) 

IF  ((MB  .EQ.  65)  .OR.  (MB  .EQ.97))  TVEN 

Al  =  KB 

A2  =  K9 

QDTG  390 
EM)IF 
IF  ((MB  .EQ.  67)  .OR.  (MB  .EQ.  99))  TVEN 

CI  =  K6 

C2  =  K7 

QOTO  390 
ELSE 

BDTD  381 
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EM)IF 
EMDIF 
IF    ((r«    .EQ.    78)    .DR.    (MB    .EQ.    110))    WEN 

GOTO  390 
ELSE 

GOrrO  380 
ENDIF 
tM)IF 

C  BY  TVE  USE  OF  TVC  CURSOR  KEYS  ON  TVE  MTERIC  KEYFVtf)  TVE  USER 

C  SELECTS  A  NEVi  SET  OF  AXIS  TO  BE  R_UI  ILU 

390  QLES  =    '     ' 

CPLL  aGTXT(28, CUES, 3, 1,335,0) 
CALL  QGTXT(28,GLJES, 3,1,320,0) 
CALL  QtJrXT(28,aJES, 3,1,305,0) 
OJES  =  (SEW  AXIS  Y  OR  N  ?' 
CALL  tX5TXT(2S, OJES, 3, 1,335,0) 
CALL  C2EEP 
CPLL  QIM<EY(I,r«) 

IF    ((r«    .ED.    89)    .OR.    (r«    .EQ.    121))    TVEN 
GOTO  400 
ELSE 

IF    ({fiB    .EQ.    78)    .OR.    (rS    .EQ.    IIO) )   TVEN 

GOTO  800 
ELSE 

GOTO  390 
EM)IF 
EM3IF 
400  QLES  =    '     ' 

CALL  QGTXT(28,aUES, 1,1, 335,0) 
CALL  aErrXT(28, QLES,  1,1,320,0) 
QLES  =    'rEW,LAST,OR  ORIGINAL    (N,L,0)' 
CALL  QGTXT(28, QLES, 3, 1,335,0) 
CALL  IA4-1-H 
CALL  QIM<EY(I,rB) 
IF    ((MB    .EQ.    79)    .OR.    (MB    .EQ.    Ill))    TFEN 
DO  410   I   =   1,N3 
XDAT(I)    =  TX(I) 
YDAT(I)    =  TY(I) 
410  CUvfTINLE 

DO  420   I    =    1,N4 
XDATl(I)    =  TXl(I) 
YDATl(I)    =  TYl(I) 
420  CONTINLE 

M  =  O 
GOTO  5 
EM)IF 

IF    ((MB    .EQ.    7<b)    .OR.    (MB    .EQ.    106))    TVEN 
XST  =   XSTl 
XFIN  =  XFINl 
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YST  =  YSTl 

'rFiN  =  YFI^a 

N5  =  M5 
N6  =  M6 
GOTD  530 
EM)IF 

IF    ((    MB    .EQ.    78)    .DR.    (MB    .EQ.    IIO) )    TVEN 
XSTl   =   XST 
YSTl   =  YST 
XFINl   =  XFIN 
YFINl   =  YFIN 
M5  =  N5 
M6  =  N6 
430  QUES  =     USE  T>t  CLFGQR  KEY  TO  MCKE' 

CPLL  aGTXT(28,ajES, 3, 1,335,0) 
QLES  =     TO  TVE  LOWER  LEFT  Hf^WD' 
CPLL.  QGrrXT(28, QLES, 3, 1,320,0) 
QLES  =     anvER,    hurt  TO  EXIT' 
CALL  QGTXT(28, QLES, 3, 1,305,0) 
CALL  QBEEP 
CALL  GLSTAT(l) 
CALL  QLTYFE(3,3) 
CALL  GU1CV(320,175) 
500  CALL  QIN<EY(I,L) 

IF    (L    .EQ.    72)    Tl-EN 
CALL  aLMD^R(0,5) 
GOTO  500 
EM)IF 

IF    (L    .EQ.    75)    THEN 
CALL  QJCIVR(-5,0) 
GOTO  500 
EM)IF 

IF    (L    .EQ.    77)    Tl-EN 
CALL  aLMD^^(5,0) 
GOTD  500 
EM)IF 

IF    (L    .EQ.    80)    THEN 
CALL  a_M»^(0,-5) 
GOTO  500 
EJVOIF 

IF    (L    .EQ.    71)    TFEN 
CALL  QLREAD(XL,YL) 
O^LL  GLM>^(320-XL,175-YL) 
CALL  QITGR(XL,YL,XST,YST) 
X2(l)    =  XST 
Y2(l)   =  YST 
CALL  QSETUP( 5,3,43,3) 
CALL  QT^«L(0,N7,X2,Y2) 
GOTO  510 
ELSE 
GOTO  430 
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EM)IF 
510  CUES  =     USE  Tl-E  CURSOR  KEY  TO  MCK/E' 

CALL  QGTXT(28,QUES,3, 1,335,0) 
QUES  =  TO  Tl-E  UPPER  RIGHT  mND' 
CPLL  QGTXT(28,QUES,3, 1,320,0) 
QUES  =  COR^eR,  HDTt  TO  EXIT' 
CPU-  QGTXT(28, QUES, 3,1,305,0) 
CPLL  QBEEP 
520  CPU.  QIM<EY(I,L) 

IF    (L    .EQ.    72)    TVCN 
CPLL  GI_MJw«(0,5) 
GOTO  520 
EM)IF 

IF    (L    .EQ.    75)    THEN 
CflLL  QUiCMR(-5,0) 
GOTO  520 
ENDIF 

IF    (L    .EQ.    77)    TVEN 
CflLL  a_ri3>^(5,0) 
GOTO  520 
EM)IF 

IF    (L    .EQ.    90)    TVEN 
CflLL  GI_MMR(0,-5) 
GOTO  520 
EM)IF 

IF    (L    .EQ.    71)    TVCN 
CPL±.  QLREPO(XR,YR) 
CPLL  Ca_rCM^(320-XR,175-YR) 
CPLL  QITGR(XR,YR,XFIN,YFIN) 
X2(l)   =  XFIN 
Y2(l)    =  YFIN 
CPiA.  aTPa_(0,N7,X2,Y2) 
CflLL  a_STAT(0) 
GOTO  530 
ELSE 
GOTO  510 
EJvDIF 
EM)IF 
530  IF    (YMPrt    .EQ.    'VOLTS')    TVEN 

GOTO  780 
E^DIF 

K6  =  1 
K7  =  N3 
DO  600   I   =  1,N3 
X   =  TX(I) 
Y=TY(I) 

IF  (((X  .GT.  XST)  -AM).  (X  .LT.  XFIN))  .AM).  ((Y  .OT.  YST) 
+        .AM).  (Y  .LT.  YFIN)))  TVEN 
K6  =  I 
GOTO  605 
ENDIF 


120 


600  CQNriMJE 

605  DO  610   I   =  K6,N3 

X   =  TX(I) 

Y  =  TY(I) 

IF  (((X  .LT.  XST)  .(jR.  (X  .ST.  XFIN) )  .OR.  ((Y  .LT.  YST) 
+       .OR.    (Y  .GT.  YFIN)))  TVEN 
K7  =  I 
GOTO  630 
ENDIF 
610        CONTINUE 
630        N5  =  K7-K6 

DO  670  I  =  K6,K7 
YDAT(I-K6+1)  =  TY(I) 
XDAT(I-K6+1)  =  TX(I) 

670        CONTINLE 
KG  =  1 
K7  =  N4 
K2  =  N4 
K3  =  1 

DO  700  I  =  1,N4 
X  =  TXl(I) 

Y  =  TYl(I) 

IF    (((X    .6T.    XST)    .PWD.    (X    .LT.    XFIN))    .PM}.    ( (Y    .GT.    YST) 
+  .«VD.    (Y    .LT.    YFIN)))    TVEN 

KS  =   I 
GOTO  715 
ENDIF 
700  CONTINUE 

715  DO  720   I   =  K8,N4 

X  =  TXl(I) 

Y  =  TYl(I) 

IF    (((X    .LT.    XST)    .OR.    (X    .GT.    XFIN))    .OR.    ( (Y    .LT.    YST) 
+  .OR.    (Y    .GT.    YFIN)))    TbEN 

K9  =   I 
GOTO  730 
EM)IF 

720  CONTINUE 

730  DO  770   I   =  K8,K9 

YDATld-KBi^l)  =  TYl(I) 
XDAT1(I-Ke+1)  =  TXl(I) 
770        CONTINUE 

N6  =  K9  -  KB 
GOTO  225 
780        N5  =  O 
N6  =  O 

DO  785  I  =  1,N3 
X  =  TX(I) 
Y=TY(I) 
IF  (((X  .GT.  XST)  .AM).  (X  .LT.  XFIN))  .AM).  ( (Y  .GT.  YST) 
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+       .fM).  (Y  .LT.  YFIN)))  TVEN 
N5  =  N5  +1 
XmT(N5)  =  TX(I) 
YDAT(N5)  =  TY(I) 
ENDIF 
7Q5        CQNTIMJE 

DO  790  I  =  1,N4 
X  =  TXl(I) 
Y=TY1(I) 

IF    (((X    .GT.    XST)    .PtO.    (X    .LT.    XFIN) )    .PM).    ( (Y    .GT.    YST) 
+  .PM).    (Y    .LT.    YFIN)))    T>€N 

N6  =  N6  +1 
XDftTl(N6)    =  TXl(I) 
YDAT1(N6)    =  TYl(I) 
EM>IF 
790  CONTIMJE 

M  =  O 
GOTO  225 
800  OJEB  =     DO  YOU  WONT  TO  PLOT' 

CPLL  QGTXT(28, OJEB, 3, 1,335,0) 
OJEB  =     THIS  GRflFH  Y  OR  N  ?' 
CfiLL  aBTXT(2B,QLJES, 3, 1,320,0) 
CPLL  C3BEEP 
CflLL  QIN<EY(I,MB) 
IF    ((MB    .EQ.   89)    .OR.    (MB    .EQ.    109))    THEN 

GOTO  830 
ELSE 

IF    ((MB    .EQ.    78)    .OR.    (MB    .EQ.    110))    TVEN 

GOTO  850 
ELSE 

GOTO  BOO 
ETDIF 
ENDIF 
830  CALL  R_aT(XDAT,YDAT,N5,XDATl,YDATl,N6,XNflrE,YNPrE, TITLE, 

+XMftX  ,  XMIN,YMftX  ,YMIN,  XST ,  XFIN,  YST ,  YFIN,  XMAJOR,YMAJOR  ,TITLE2, 
-tFACT,NflrE) 

IF    (YNPrt    .EQ.    'V^1.TS' )    TFEN 

GOTO  390 
EM)[F 
IF    (YNflhE    .EQ.        Dv/Di)    TVEN 

GOTO  380 
E^OIF 
850  DO  860   I   =  1,2000 

TX(I)=0 
TY(I)=0 
TXl(I)    =  O 
TYl(I)    =0 
XDAT(I)    =  O 
XDATl(I)    =  O 
YDAT(I)    =  O 
YDATl(I)    =  O 
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860  CONTINLE 

CflLL  QEMjDE(16) 
CPU.  rSCiDE(3) 

900  FDRMAKAl) 
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THIS  SLSFOJTIhE  RJDTS  THE  GRf^FHS  AS  SHDMM  ON    fVE  SCREEN  DURING 
TVC  CPLLING  SLff«XmrsE  GFl^FHl. 

THE  MAJORITY  OF  TVE  PRDGRfW   IS   IN^/DLSO)   IN  SETTING  LF  flRGLJE^ENTS 
TO  BE  PASSED  TO  SLS^DLfTir^  CO^^'AI^eD  WITHIN  T>E  R_OTMATICS 
LIBRWY.    THIS  MANUAL  SHDLLD  BE  CONSULTED   IN  ORDER  TO  DETEPMINE 
THE  USE  OF    rVE  PR3UE^ENTS  OTVER  1V¥^  THOSE  CEFir^ED  BELOW 

XDAT      :    THE  ARRAY  CONTAINING  T>E   X  V¥^UES  OF  TVE  FIRST  CURVE 

TO  BE  PLOTTED. 
YDAT      :    Tl-E  FIRST  CURVE'S  Y  VPLUES 
XDATl    :    TVE  ARFtf^Y  CDNATININB  IVE.  Y  VALUES  OF  TVE  SECOND  CURVE 

TO  BE  PLOTTED. 
YDATl    :    TVE  SECOND  CURVE'S  Y  VVIjUES 
TITLE    :    TVE  TITLE  OF  THE  PLOT 
"    TLE2:    AN  DISCRIBTIDN  OF  TTE  GRPPH 
XNP^E    :    TVE   X  AXIS  TITLE 
SATf^E:    TVE  Y  AXIS  TITLE  PROPERLY  POSITION  PND  CATENTATED 

SUBRDLm^E  PLDT(XDAT,YDAT,^FT1, XDATl, YDATl, ^PT2,XNA^E,YNA^E 
+  , TITLE, XMAX,XMIN,YMAX,YMIN,XST,XFIN,YST,YFIN,XliAJOR,YMAJOR 
+  ,TITLE2,FACT,    NATE) 

REAL*4   XDAT ( 2000), YDAT (2000), XDATl (2000), YDATl (2000) 
INTEGER*2  JCOL(IO)  ,JROW(  lO)  ,M^1,NFT2 
INTEGER   13, 14, 15, 16, 17, 18, 111, 112, re 

CHARACTER*8  NAI-E,YNPIt«5,XNfil-E«ll,TITLE«15,SAr'M£«12,TITLE2*32 
+,    FACT«6,A4*1,    QUES«28 
A4  =    '     • 

SATFLE  =  YNATE/ZAA/ZFACT 

XORG  =   XST 

YORG  =  YST 

DTX   =  ABS(XMAX-XMIN) 

DX   =  DTX    /MAX(^PT1,^FT2) 

DTY  =  ABS(YMAX-YMIN) 

lOPT  =  O 

M<MT  =   5 

NDOTS  =  O 

LABEL  =  -1 

ivoa:  =  3 

ISIZE  =  1 

ISTEP  =  2 

IPORT  =   1 

CALL   ZINIT(  ISIZE,  ISTEP,  IPORT, BATCH) 

JCL  =  20 

JCR  =  620 

ICL  =   10 

ICR  =  320 

13  =   12 

14  =   150 
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15  =  236 

16  =  336 

17  =  340 

18  =  322 

111  ^  386 

112  =  182 

CALL  ZC0N^(JCL,ICL,16) 
CALL  ZCa^(JCR,ICR,16) 
CALL  ZCDNV(  13, 14,16) 
CALL  ZC0NV(I5,I6,16) 
CALL  ZCGNV(  17, 18,16) 
CPU-   ZCONV( 111,112,16) 
XMINl  =XMIN  -.2*DTX 
XMRXl  =  XMAX+   -lt(DTX) 
YMAXl  =  VT1AX  +  .ItDTY 
YMINl  =  YMIN  -  .1«DTY 
MDGE  =  2 

CPLL  ZR.OT(JCL,JCR,ICL,ICR, XMINl, XMAXl, YMINl, Y?1AX1,XCF!G,YCRG) 
CALL  ZFfTOI(XST,YST,Ll,L2) 
CALL  ZRTOI(XFIN,YFIN,L3,L4) 
CALL  ZSETUP(-1, 1,0,1) 
MINOR  =  1 

CALL  ZXAXIS(XST,XFIN,XMP^CR,MINCF?, LABEL, M3EC) 
CALL  ZYAXIS(Y5T,YFIN,YMAJCR, MINOR, LABEL,NDEC) 
CALL  ZXTICS(lO,JC0L,JRDW) 
CALL  ZYTICS(10,JCCL,JRDW) 
CALL  ZLI^E(L1,IJ2,L1,L4,1) 
CALL  ZLI^E(L1,L4,L3,L4,1) 
CALL  ZLir^C(L3,L4,L3,L2,l) 
CALL  ZLIIVE(L3,L2,L1,L2,1) 
CPLL   ZRJ 
QUES  =  '  ' 

CALL  aGTXT(28, QUES, 3, 1,335,0) 
CALL  GGTXT(2B, QUES, 3, 1,320,0) 
CALL  QETXT(28, QUES, 3, 1,305,0) 
OJES  =  '  IF  DESIFtD  INSERT  A  bEM   PEN' 
CALL  aGTXT(28, QUES, 3, 1,335,0) 
QUES  =  Tl-EN  PFE35   ANY  KEY- 
CALL  QGTXT(2B, CUES, 3,1,320,0) 
CALL  QBEEP 
CALL  QIN<EY(I,M8) 
CALL  ZTABL(1,NPT1,XDAT,YDAT) 
CALL  ZRJ 

OJES  =  '  IF  DESIRED  INSERT  A  f^W  PEN' 
CALL  Gm"XT(28,aiJES, 3, 1,335,0) 
OJES  =  TVEN  PRESS  ANY  KEY- 
CALL  QGTXT(28,aJES,3,l,320,0) 
CALL  QBEEP 
CALL  QIM<EY(I,r«) 
CALL  ZTABL(1,^FT2,XDAT1,YDAT1) 
CALL  ZPU 
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CflLL  QIM<EY(I,I-B) 

CUES  =    '  IF  DESIRED   irGERT  A  iVEW  PEN' 

CPLL.  QGTXT(28,GLE5, 3, 1,335,0) 

OLES  =      TVEN  PRESS  PNY  KEY' 

C^LL  C3GTXT(2a,aLJEB, 3, 1,320,0) 

CPU-  QBEEP 

CPLL   ZPTXTA(ll,XNflre,l) 

CPLL   ZSFrLP(2, 1,0,1) 

CfiLL  ZIW(  1,1, 1760,1400,0) 

CALL  ZPTXT(15,TITLE,1,I5,I6) 

CPLL  ZPTXT{8,Nflrt,  1,111, 16) 

CALL   ZPTXT(32,TI"n_E2,l,I5,I8) 

CALL   ZEXTXT(0,0,0,0,0,1) 

CALL   ZPU 

CALL   ZPA(I3,I4) 

CALL   ZLABELC 12, SAMPLE) 

CALL   ZFINIS 

END 
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THIS   IS  Tl-E  FIRST  SLBROLfTirC  CALLED  W-EN  USING   THE  CUBIC 
SPLI^E  rtTHTO. 

MDATA   :    TVE  PfV^Y  CDhfTAININB  POTENTIPL  PM)  LOB  CLFiFENT  ESVEITY 
XEAT      :    TVE  PF¥ifi>W  CG^ITAINI^G  TVE  LOG  CURFCNT  V-V^LUES  FDR  TVE 

CATHDDIC  BRfiNCH 
YDAT      :    ThE  P«RAY  CC^r^AINI^G  TVE  DERIVATIVES  FOR  TT-E  CATHDDIC 

BRANCH 
XDATl    :    S<^«t  A  XDAT  BUT  FDR  T>E  AM3DIC  BRANCH 
YDATi    :    SATE  AS  YDAT  BUT     FOR  TVE  PNODIC  BRANCH 
DIV        :    TVE  ARFtf^Y  CCNTAINING  TVE  DERIVATIVES  UHICH   IS  BEING 

RETURN  FROM  SUBSEQUENT  CALLED  SUBRDUTI^ES 

TVE   INTEGER  VARIABLES  ACTUALLY  USED   IN  THIS  SLBRDUTIIVE  ARE 
USED  TO  DEFINE  TVE  REGIOTS  WITHIN  BDTTH  BR^1M>ES  AVAILAtLE 
FOR  USE 


SUBRDLrriNE  SLCFEl(riDATA,N,J,NfirE,Al,A2,Cl,C2) 

REAL*4  XDAT ( 200O ), YDAT ( 2000), XDATl (2000), YDATI (2000) 

FEAL»8  MDATAC 2000,2)  ,DIV(1000) 

INTE6ER<2   I  ,J,N,K,MI  ,N3,N4,N1,K1,A1,A2,C1,C2 

CHARACTER  NAf«E:«8,YNArE»5,TITLE2«32,TITLE3«50 

Nl   =  N 

YNfirt  =     Dv/Di' 

TITLE2  =     TAFEL  SLOPES:   CUBIC  SH_I^E• 

TITLE3='     • 

C     TVE  CATHDDIC  BRPNCH  IS  Cl-ECKED  FDR  AN  INCREASE  IN  LOG 
C     CURRE^fT  DENSITY  WHICH  COULD  CAUSE  TVE  CALCULATION  OF  A 
C     DERIVATIVE  OF  IM^  INTE  SLOPE 

DO  20  I  =  J-1,2,-1 

IF  (M>ATA(I-1,2)  .LE.  M)ATA(I,2))  GOTO  25 
20    CONTIMJE 

K  =  MAX((J-250),1) 
GOTO  26 

25  K  =  I 

26  MI  =  J-9 
Jl  =  K 

N3  =  ^«S(K-MI) 
M  =  7 

CALL  CSR_IN(MDATA,N,M,DIV,K,N3) 
DO  100  I  =  l,N3-3 
XDAT(I)  =  rDATA(J-4-I,2) 
YDAT(I)  =  DIV(N3-I+1) 
100   COMTINUE 
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C  THE  fiMDDIC  BFtf^fCH    IS  CHECKED  FDR  LUG  OJRf^NT  WfiUJES  WHICH 

C  GjULD  cause  TFE  CAUaJLATION  CF  A  CGRIVATIV^  GF    IM^'INITE  SLjDPE 

DO  120   I    =  J+l,r41 

IF    (r«ATO(I,2)    -LE.    ri>ATA(  1-1,2))    GOTO  125 
120        COvn'IMJE 

K  =  MIN((J+250),(N1)) 
GOTO   126 

125  K  =   I 

126  Jl   =  J 

m   =  K  -  9 
N4  =  ABS(MI-J) 

CALL  C9PLIN(MDATA,N,ri,DIV,J,N4) 
DO  200   I   =   1,N4 
XDATl(I)    =  MDATA(Jl+I+4,2) 
YDATl(I)    =   DIV(I) 
200        CCTiTINLE 

D^LL  GF¥¥^l(XDAT,YDAT,N3,J,IMflrt,XDATl,YDATl,N4,YNfirt,TITl_E2, 
+A1,A2,C1,C2,TITL£3) 
E^© 
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THIS  SUff«Xm^E  sets  LF  D-E  MATTRIX  CF  COEFFICIQVfrS  TD 
BE  SOLVED   IN  ORDER  TD  DETERMI^E  71-E  DERIVATIVE  USI^G  THE 
CLSIC  SR_irvE  rtTTHOD, 

MDATA    :    ThE  Pffif^Y  CC^fTAINI^G  POTENTIAL  AM)  LOG  CURRENT  DE^eITY 
CrV^T      :    T>E  7x7  MATTRIX  CINTAININB  TVE  COEFFICIENTS 
YMAT      :    TFE  7x1   VECTOR  CONTAINING  TVE  RIGHT  HflND  SIDE  OF 
TFE  SYSTEM  CF  EGUATIOTvE.    W-EN  RETLRSED  FROM  TFE 
SOLVING  SLBROUTirc:   IT  CONTAINS  TFE  SOLUTION. 

SUBRGLJTIhE  CSR_INE(M)ATA,J,M,DIV,K,MI) 
HEPL«8   XDAT(IOOO),   YDAT(IOOO)  ,M)ATA(2000,2)  ,YMAT(7) 
REAL*8  CMAT(7,7),AI,BI,CI 
REPL*8  DIV(IOOO) 
irrTEGER*2  N,MI,K,J 
DO   lO   I    =   1,7 
DO  5  J  =   1,7 

CMATd.J)    =  O.DO 
5  CCNTIMJE 

lO  CONTINUE 

M  =  7 
K  =  K-1 

DO  200  N  =   1,MI 
20  K  =  K  +   1 

CMAT(1,1)    =  2.D0«(rf)ATA(K+2,2)-rCATA(K,2)) 
CMAT(1,2)   =  M)ATA(K+2,2)    -  MDATA(K+1,2) 
CMAT(7,7)    =  2.D0*(M)ATA(K+e,2)-W)ATA(K+6,2)) 
CMAT(7,6)    =  MDATA(K+8,2)-M3ATA(K+7,2) 

YMAT(1)=  ((M)ATA(K+2,l)-M^TA(K+l,l))/(n3ATA(K+2,2)-MDATA 

+  (K+1,2)))    -    ((fDATA(K+l,l)-M)ATA(K,l))/(M)ATA(K+l,2)-M)ATA 
+(K,2))) 

VMAT(l)    =  YMAT(1)*6. 

YMAT(7)=         ((MDATA(K+8,l)-MDATA(K+7,l))/(M)ATA(K+8,2)-f-CATA 
>(K+7,2)))    -    ((^DATA(K+7,l)-^OATA(K-^6,l))/(MDATA(K+7,2)-^CATA 
+(K+6,2))) 
YMAT(7)    =  YMAT(7)«6. 
DO   100   I   =  2,6 
K  =  K  +1 

CMAT(I,I)  =  2.DO«((M}ATA(K+2,2)-fDATA(K,2))) 
CMAT(  1,1+1)  =  M>ATA(K+2,2)-MDATA(K+i,2) 
CMAT(  1,1-1)  =  M)ATA(K+1,2)  -  M)ATA(K,2) 
YMAT(I)  =  ((M>ATA(K+2,1)-M)ATA(K+1,1))/CMAT(I,I+1)) 
-»-( (M>ATA(K+1,1  )-M)ATA(K,l )  )/CMAT(  1,1-1)) 
YMAT(I)  =  VMAT(I)«6. 
100   CCNTIMJE 
K  =  K  -  5 

CALL  LIN5Y1(CMAT,YMAT,M) 
AI  =  YMAT(5)  -  YMAT(4) 
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AI  =  AI/((6.D0)*((MDATA(K+4,2)-f'OATA(K+3,2)))) 
Bf  =  \T1AT(4)/2.DO 

CI  =  (r'5)ATA(K+4a)-MDATA(K+3,l))/(rT)ATA(K+4,2)-riDftTA(K+3,2)) 
CI  =  CI-(l/6)«(((2<YMAT(4))+YMAT(5))*(MDATA(K+4,2)-f1DATA(K+ 
+3,2))) 
DIV(N)  =  CI 
200   CONTINUE 
END 
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THIS  SUBROUTIJ^  SOLVES  Tl-E  7x7  SYSTEM  fJF  EOL^TICrMS  GBVEFATED 
BY  rJEPLIN.FOR.      TVE  FRDGR^  WftS  CFIGINPLLY  FRCM  THE  TEXTBOOK 
"FCFTfWN  77   for  Enqineers  and  Scientists"   BY  NVHDFF  ff€> 
LEESTMA,    PAGES  264-265.    ThE  MAIN  MCDIFICATION   IS  TFE  ELIMIf^TICJvJ 
OF  TFE  VPRIABLE  DirGvJSICNING  OF  ARRAYS  PRESENT    IN  TVE  ORIGINAL 
PROGRAM. 

SLJBRDL^■I^E  LIN5Y1(CCMAT,YMAT,N) 
INTELER  N.PIVCrr 

REPL*8  ALX3(7,8),    CDMAT(7,7),    YMAT(7)  ,X(7)  ,MJLT,TE]^ 
DO  '20   I    =   1,   N 
DO   10  J   =   1,N 

ALG(I,J)    =  CCriAT(I,J) 
10  CONTIMJE 

20  COrfTIMJE 

DO  30   I   =  1,N 

ALJG(I,N*-1)    =  YMAT(I) 
30  CONTIMJE 

DO  70   I   =   1,N-1 

IF    (ALG(I,I)    .EQ.    O.)    Tl-EN 
PIVOT  =  O 
J   =   I+l 
35  IF    ((PIVOT    .EQ.    O)    -AND.    (J    .LE.   N) )    Tl-EN 

IF    (ALJG(J,I)    .NE.    O. )    PIVOT  ^  J 
J   =  J+l 
GOTO  35 
ENDIF 

IF    (PIVOT    .EQ.    O)    TVEN 
PRINT*,    'SI^GULAR• 
STOP 
ELSE 
DO  40  J   =   1,N+1 
TETf^  =  ALG(I,J) 
ALJG(I,J)    =  ALJG(PIVDT,J) 
ALJ6(PIVDT,J)    =  TEl^ 
40  CONTINUE 

ENDIF 
ENDIF 
DO  60  J   =    I+1,N 

MULT  =  -ALJG(J,I)/ALJG(I,I) 
DO  50  K  =   I,N*-1 

ALJG(J,K)    =  ALJG(J,K)    +  riJLT*ALG(  I  ,K) 
50  CONTINUE 

60  CONTINUE 

70  CONTINUE 

X(N)    =  AUB(N,N+1)/AUG(N,N) 
DO  90  J   =  N-1,1,-1 
X(J)    =  AL6(J,N<-1) 
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DO  80  K  =  J>-1,N 

X(J)    =  X(J)-(^^UG(J,K)«X(K)) 
QO  CXlSfriMJE 

X(J)    =  X(J)/«JG(J,J) 
90  CON^riMJE 

DO  lOO    I   =  1,N 
VriATd)    =   X(I) 
100  aJMTIMJE 

EM) 
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«  THIS  SUBRDUTIIVE  CXJNJERTS  A  MJhERIC  V^^IBLE    (X6)    OF  FCH^T 

«  +/-  XX.XXXX   TO  A  CHflRACTER  VflRIBLE  WHOSE  WIDTH   IS  8  CHPRACTEFS. 

«  THIS   IS  ^ECESSARY  AS  Tl-E  PLDTMATICS  LIBRflRY    IS  NOT  EQUIPPED  WITH 

*  A  SUBROUTINE  FOR  R-OTTINE  MJTtRIC  VARIBLES. 


SUBROUTirvE  CaMV(A7,X6) 
irvnEGER  01,02,03,04,05,06 

C1,C2,C3,C4,C5,C:6,C7,C8,A7»8 


CHflRAUItK*!  A6(0:9), 

C1,C2,C 

A6(0)  =  0' 

A6(i)  =  1' 

P^(2)  =  '2' 

A6(3)  =  '3' 

A6(4)  =  '4' 

A6(5)  =  '5' 

A6(6)  =  '6* 

A6(7)  =  '7' 

A6(8)  =  8' 

A6(9)  =  '9' 

CI  =  '  ' 

C2  =  '  ' 

C3  =  '  ' 

C4  =  •  ' 

C5  =  '  ' 

C6  =  •  ' 

C7  =  '  ' 

IF  (X6  .LT.  0.0)  CI  = 

=  '-• 

X6  =  ARR(X6) 

XI  =  X6/10. 

01  =  IFIX(Xl) 

X2  =  X6-(10.»01) 

02  =  IFIX(X2) 

X3  =  X6  -  (lO.WDl)  - 

02 

X3  =  X3/.1 

03  =  IFIX(X3) 

X4  =  (X3-a3)/.l  +  .OOl 

04  =  IFIX(X4) 

X5  =  (X4-04)/.l 

05  =  IFIX(X5) 

X7  =  (X5-05)/.l 

06.  =  IFIX(X7) 

DO  10  I  =  0,9 

IF  (01  .ED.  I)  C2 

=  A6(I) 

IF  (02  .ED.  I)  C3 

=  A6(I) 

IF  (U3  .ED.  I)  C4 

=  A6(I) 

IF  (04  .ED.  I)  C5 

=  A6(I) 

IF  (05  .ED.  I)  Ch 

=  A6(I) 
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IF    (06    .ED.    I)    C7'  =  A6(I) 
10  COrvfrKVLE 

ca  =  '.' 

A7  =  C1//C2//C3//C8//C4//C5//C6//C7 

RETURN 

E^JD 
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«         THIS  ^^FO-^"I^E  sets  sort  of  tj-e  a^mfiCJUi  strimbs  used 

»  W-EN  USI^G  Tl-E  CENTTVi-  DIFFEI^NCE  ItlT-DD 

SL»«JTI^E  3JCFE(n)ATA,N,J,N«*E,Al,A2,Cl,C2,Nl,N2,N3,N4) 
REPL«4  XDAT(2000),YDAT( 2000 ),XDAT1( 2000 ),YDATl(2O0O) 
REAL»8  MDATA(2000,2) 
INTEGERrZ  I,J,N,A1,A2,C1,C2 

CHflFWCTER  ^W^t«8,YNfirE*5,TITLEZ•32,    LIST«50,TITLE3«50 
YNPME  =    Dv/Di' 

TIT1_E2  =     TAFEL  SLOPES:    CENTPAL  DIFFEFOCE' 
TITLJE3='     ' 

CALL  DATADEL(MDATA,J,YDAT,YDAT1,XDAT,XDAT1,M2,N1) 
CALL  ffV:FHl(XDAT,YDAT,N3,J,rv*2rE,XDATl,YDAri,N4,YIS¥:»t,TITLE2, 
+A1 ,A2,C1 ,C2,TITLE3) 
END 
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THIS  SLe^OL^"I^E  ccnducts  mj-ekicp*.  diffepentiation  lbing 

THE  4  POINT  CENTTV^   DIFFERNCE  fETHDD. 

r«o=>TA    =  p«RAY  a3vrrAINI^e  pdtentipl  and  log  i  reiveiTY 

SLCFE      =  PfT^AY  CCNTAINING  De/D(LDGI )    ATC  LOG    I 

FSLDFt  =  CENTRAL  DIFFEPTCE  DIFFEFENTIATION 

ECORRl   =  LPPEHLIMIT  CF  ATVCDIC/CATHDDIC  DISTRIBLTTIGN 

EC0FR2  =  LC]l«EFLIMIT  OF  PMX>IC/CATHDDIC  DISrmiBJTION 

N  =  MJMBER  OF  13ATA  POINTS 

J  =  LOCATION  WITHIN  fDATA  OF  ECORR 

I  =  CAT>a)IC  BR=M>I  LOOP  COUNTER 
K               =  ftfJDIC  BFV^NCH  LOOP  COUNTER 
H  =  DIFFERENTIAL  MPLUE   IN  POTENTIAL 

II  =  VV^LE  OF  LOG    I    AT  E-2H 

12  =  VALUE  OF  LOG    I    AT  E-fl 

13  =  V«JJE  OF  LOG   I   AT  E 

14  =  V/flLUE  OF  LOG  I  AT  E+H 

15  =  VPLUE  OF  LOG  I  AT  E+2H 
DEDI  =  SLOPE  VALUE  OF  13 

SUBRDLrrilVE  DATADEL  (W)ATA,J,YDAT,YDATI,XDAT,XDAT1,N2,N1) 

REAL«a  M3ATA(2OOO,2),Xl,X2,X3,X4,X5,H,FSL0PE,DEDI 

F=EAL«4,  XDAT1(1000),YDAT1(1000),  YEAT(IOOO)  ,XDAr(1000) 

INTEGER«2  J,K,I 

FSLOPE  (H,Xl,X2,X4,X5)=(Xl+(a.IX)«X4)-(8.D0«X2)-X5)/(12.D0«H) 

C      PERFORMING  CALCULCATIONS  ON  CATHODIC  BRANCH 

K  =  O 

DO  100  I  =  J-2,Nl+2,-l 

K  =  K+1 

H  =  rT)ATA(I,l)  -  MDATA(I-1,1) 

XI  =  rCATAC 1-2,2) 

X2  =  MDATA( 1-1,2) 

X3  =  riDATA(I,2) 

X4  =  MDATA( 1+1,2) 

X5  =  MDATA( 1+2,2) 

DEDI  =  l.DO/FSL0FE(H,Xl,X2,X4,X5) 

YDAT(K)  =  DEDI 

XDAT(K)  =  riDATA(I,2) 
lOO    CONTINUE 
1=0 

C  FEFFDR1ING  CALCULATIONS  ON  TT-E  ANODIC  BRflNCH 

200  DO  300  K  =  J+2,    N2-2 

1    =    1+1 
H      =  riDATA(K+l,l)    -  M}ATA(K,1) 
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XI  =  MDATA(K-2,2) 

X2  =  nDATA(K-l,2) 

X3  =  MmTA(K,2) 

X4  =  riDATA(K+l,2) 

X5  =  MDATA(K+2,2) 

DEDI  =  l.DO/F5LXPE(H,Xl,X2,X4,X5) 

YDATl(I)  =  DEDI 

XDATl(I)  =  f-OATA(K,2) 
300    CDNTINUE 
400    RETURN 
EIvO 
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